THE 


AERONAUTICAL 
QUARTERLY 


A Journal devoted to Aeronautics and the Allied Sciences 


Volume XI NOVEMBER 1960 Part 4 


CONTENTS 
Graphite: A Future Structural Material A. J. Kennedy 


Boundary Layers over Infinite Yawed 
Wings J. C. Cooke 


A Small-Bore Pick-Up for the Measure- R. E. Franklin and 
ment of Fluctuating Pressures R. B. Archbold 


The Aerodynamic Forces Associated 
with Harmonic Motion of Slender 
Wings on Stationary Bodies of 
Revolution R. D. Milne 


The Reliability of Multiple Rocket 
Engines Gillian Ruston 


The Quasi-Cylinder of Specified Thick- 

ness and Shell Loading in Supersonic 

Flow H. Portnoy 
Review of ‘Progress in Solid 

Mechanics,” Vol. I, edited by I. N. 

Sneddon and R. Hill W. S. Hemp 


Index to Vol. XI 


LONDON 


ROYAL AERONAUTICAL SOCIETY 
4 HAMILTON PLACE WI 


Ja 


THE AERONAUTICAL QUARTERLY 


Editorial Board 


Dr. A. M. BALLANTYNE, T.D., B.Sc., Ph.D., F.R.Ae.S., Hon.F.C.A.I. 
F.LA.S. (Chairman) 


Mr. HAYNE CONSTANT, C.B., C.B.E., M.A., F.R.S., F.R.Ae.S. 
Sir HAROLD ROXBEE COX, D.Sc., D.LC., Ph.D., F.R.Ae.S., F.LA.S. 
Sir WILLIAM FARREN, C.B., M.B.E., M.A., D.Sc., F.R.S. 


Hon.F.R.Ae.S., Hon.F.LA.S. 


Sir ARNOLD HALL, M.A., F.R.S., F.R.Ae.S. 
Sir BENNETT MELVILL JONES, C.B.E., A.F.C., M.A., F.RS., 


Hon.F.R.Ae.S., Hon.F.LA.S. 


Sir ALFRED PUGSLEY, O.B.E., D.Sc., F.R.S., F.R.Ae.S. 
Sir GEOFFREY TAYLOR, M.A., D.Sc., F.R.S., Hon.F.R.Ae.S. 
Dr. P. B. WALKER, C.B.E., M.A., F.R.Ae.S. 


Aerodynamics 


(a) General 

Sir Leonarp Bairstow 
Mr. D. Kertn-Lucas 
Prof. W. A. Mair 

Mr. L. F. NICHOLSON 
Dr. R. C. PANKHURST 


Dr. B. THwaltes 


(b) Fluid Motion 
Dr. G. K. BaTCHELOR 


Prof. J. H. 
Prof. H. B. Squire 
Prof. G. TemMPLe 


(c) Stability and Control 


Mr. L. W. Bryant 

Mr. S. B. Gates 

Mr. H. B. Invino 

Mr. G. H. Leg 

Mr. M. B. Morgan 

Mr. H. H. B. M. THomas 


Helicopters and 
Propellers 


Prof. J. A. J. Bennett 
Mr. P. Bagrr 

Mr. L. G. Fatrnurst 
Mr. R. Harner 

Dr. G. S. Histor 
Capt. R. N. 


Editor—E. C. Pike, M.A., A.F.R.Ae.S. 


Editorial Referees 


Instruments and 


Electrical Equipment 
Dr. G. E. Bairsto 

Sir Georce GARDNER 

Dr. C. S. HupDsON 

Mr. C, N. Jaques 

Mr. F. W. MEREDITH 


Aircraft Design 
. R. 
. R. M. CLARKSON 
. A. N. CLIFTON 
. H. H. Garpner 
. H. KNOWLER 
. E. W. Perre 


Dr. A. E. RUSSELL 
Mr. R. S. STAFFORD 
Mr. W. H. STEPHENS 


Radio and Radar 

Mr. J. F. ATHERTON 

Mr. B. G. Gates 

Dr. R. L. Rose 
Dr. A. M. UTTLey 

Sir Roperr Watson-Watt 


Aviation Medicine 


Dr. K. G. 
Dr. A. BUCHANAN Barspour 


Meteorology 

Sir Davip Brunt 
Mr. J. K. Harpy 
Prof. P. A. SHEPPARD 


Sound 

Mr. W. MAKINSON 
Mr. F. B. Greatrex 
Mr. N. FLEMING 


Aircraft Propulsion 
Mr. A. D. BAXTER 
Mr. A. V. CLEAVER 
Sir A. H. Roy FEDDEN 


Mr 
Prof. O. A. SAUNDERS 
Prof. A. G. SMITH 


Air Transport 


. A. C, CAMPBELL Orne 
. P. G. MASEFIELD 

. N. E. Rowe 

. B. S. SHENSTONE 

. W. Tye 


Aircraft Accessories 


Mr. R. GraHaM 
Mr. C. D. HOLLAND 


Performance 


Flight Testing 


Mr. R. R. Duppy 
Mr. P. A. HUFTON 


Structures 


Prof. J. H. Arcyris 
Mr. A. J. Barretr 

Mr. A. C. Brown 

Mr. H. L. Cox 

Mr. D. J. Farrar 
Prof. W. S. Hemp 
Prof. R. Hitt 

Mr. H. B. Howarp 
Col. H. T. Jessop 

Mr. E. D. KEEN 

Mr, F. A. Kerry 

Dr. D. M. A. LeaGetr 
Mr. R. L. Licxigy 
Prof. A. J. SUTTON PipPaRD 
Prof. S. C. REDSHAW 
Mr. R. H. SANDIFER 
Sir Richarp SOUTHWELL 


Fuels and Oils 


Mr. E. L. Bass 

Dr. J. W. DrinkwaTer 
Prof. F. H. GARNER 
Prof. W. E. Garner 
Dr. K. E. W. RIDLER 
Dr. C. G. WILLIAMS 


Materials 


. L, ROTHERHAM 
Dr. H. Sutron 
Dr. C. Sykes 


Vibration and Flutter 
Mr. E. G. BrOaDBENT 
Prof. A. R. COLLAR 

Prof. W. J. DUNCAN 

Mr. H. TEMPLETON 

Mr. M. O. W. WOLFB 


Dr. S. Hooker 
Mr. A. R. HOWELL 
Mr. W. H. LINDSEY 
Mr. P. Lioyp 
Mr. E. F. Ree Dr. A. W. Morey 
Prof. E. J. RICHARDS Dr. E. S. MOoULT 
Mr. N. S. Muir 
Mr. F. Nrxon 
Mr. F. M. Owner 
Dr. D. W. HOLpsr Mr. F. H. PoLiicutr 
Prof. L. HowarTtH 
Mr 
Mr 
Mr 
Prof. G. N. WarD 
Prof. A. D. Youno 
Mr. H. G. Conway PF 
Str Dowry Dr. L. 
Dr. Aum 
Prof. N. F. Morr 
Mr. D. R. NEWMAN 
Mr. A. H. Yates 
Mr. E. T. JONES 
Mr. S. Scott 
Mr. H. L. STEVENS 


Graphite: A Future Structural Material 


A. J. KENNEDY, Ph.D., A.M.LE.E., F.1.M., F.Inst.P. 


(The College of Aeronautics, Cranfield) 


SUMMARY: Graphite is a material with unique high-temperature properties. Its 
strength and modulus increase with temperature up to about 2000-2500 C., in which 
range it exhibits a higher strength than any known material. The development ot 
new types of graphite (such as pyrolytic graphite, deposited from the vapour) and 
impregnated graphites, are both significant to future technology. The highly- 
oriented pyrolytic graphites in particular, with their marked thermal and mechanical 
anisotropy, and their superior mechanical properties (60 ~ 103 Ib./in.? strength at 
2750 C., with 60 per cent ductility) offer distinct possibilities, particularly in 
composite structures. Graphite sublimes directly from the solid to vapour phase at 
3700 C. at atmospheric pressure. As the sublimation energy is high, it is, at least in 
theory, a very efficient ablating material. It also has a high creep resistance. Sur- 
face reactions limit its value, particularly in oxidising environments, and coatings, 
or other methods of surface protection, are of high importance. The modulus of 
graphite is low, being about one-half of that for lead at room temperature, 
increasing with temperature to a maximum in the region of 2200 C. The properties 
vary in detail from one type of graphite to another, and with the direction in which 
measurements are made. Graphite has a low density (about 2:0) and is com- 
paratively cheap, so that its low temperature limitations (zero ductility at room 
temperature) may be offset, at least in part, by a re-appraisal of the design factors. 
It can be joined and fabricated fairly readily, and appears to be a material of great 
value to advanced high-temperature projects in the future, where composite 
assemblies of various kinds will be essential. 


1. Introduction 


Until very recent times, temperature questions were not of any very serious 
concern in aircraft structures. Materials for structural applications were evolved 
largely on the basis of strength weight ratios and stiffness, although higher speeds, and 
increased dynamic stressing, brought fatigue resistance into such prominence that 
this factor has exercised, in latter years, a decisive influence on both design and 
materials utilisation. While fatigue has also been of concern in the development of 
propulsion units, deformation of engine components by creep at the operating 
temperatures has been the controlling materials factor, and the need for materials 
capable of sustaining working stresses at temperatures around 1000 C. has provided a 
considerable incentive to recent metallurgical research and development in this 
direction. 
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Graphite: A Future Structural Material 


A. J. KENNEDY, Ph.D., A.M.LE.E., F.I.M., F.Inst.P. 
(The College of Aeronautics, Cranfield) 


SumMARY: Graphite is a material with unique high-temperature properties. Its 
strength and modulus increase with temperature up to about 2000-2500°C., in which 
range it exhibits a higher strength than any known material. The development of 
new types of graphite (such as pyrolytic graphite, deposited from the vapour) and 
impregnated graphites, are both significant to future technology. The highly- 
oriented pyrolytic graphites in particular, with their marked thermal and mechanical 
anisotropy, and their superior mechanical properties (60 x 103 Ib./in.2 strength at 
2750°C., with 60 per cent ductility) offer distinct possibilities, particularly in 
composite structures. Graphite sublimes directly from the solid to vapour phase at 
3700°C. at atmospheric pressure. As the sublimation energy is high, it is, at least in 
theory, a very efficient ablating material. It also has a high creep resistance. Sur- 
face reactions limit its value, particularly in oxidising environments, and coatings, 
or other methods of surface protection, are of high importance. The modulus of 
graphite is low, being about one-half of that for lead at room temperature, 
increasing with temperature to a maximum in the region of 2200°C. The properties 
vary in detail from one type of graphite to another, and with the direction in which 
measurements are made. Graphite has a low density (about 2-0) and is com- 
paratively cheap, so that its low temperature limitations (zero ductility at room 
temperature) may be offset, at least in part, by a re-appraisal of the design factors. 
It can be joined and fabricated fairly readily, and appears to be a material of great 
value to advanced high-temperature projects in the future, where composite 
assemblies of various kinds will be essential. 


1. Introduction 


Until very recent times, temperature questions were not of any very serious 
concern in aircraft structures. Materials for structural applications were evolved 
largely on the basis of strength/weight ratios and stiffness, although higher speeds, and 
increased dynamic stressing, brought fatigue resistance into such prominence that 
this factor has exercised, in latter years, a decisive influence on both design and 
materials utilisation. While fatigue has also been of concern in the development of 
propulsion units, deformation of engine components by creep at the operating 
temperatures has been the controlling materials factor, and the need for materials 
capable of sustaining working stresses at temperatures around 1000°C. has provided a 
considerable incentive to recent metallurgical research and development in this 
direction. 
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The construction of aircraft capable of speeds in the region of Mach 2-3 raises a 
whole series of new issues as far as materials are concerned, the most obvious being 
the achievement of skin temperatures of about 100-150°C. This is precisely the 
temperature range in which the strengths of the conventional aluminium alloys fall off. 
Creep also starts to become a significant consideration at these temperatures. It may 
well be possible, by the best use of materials, together with suitable cooling techniques, 
to retain aluminium alloys as the principal structural materials for such aircraft, but 
there is no doubt that future technology must be prepared to embark upon quite new 
ventures in materials if advanced machines, aeronautical or astronautical, are to be 
successfully brought into being. Were it not for the temperature factor, it is highly 
likely that the present range of available materials would be adequate, even though 
improvements are naturally always desirable. The existence of high temperatures, 
and of rapid rates of heating, transforms the whole scene, not merely because 
temperature-resistant materials are now required where conventional alloys once 
sufficed, but because of a whole range of secondary issues which involve weight, 
compatibility, joining, deformation, conductivity and expansion considerations, as 
well as many others. 

Examination of the range of possible materials, including those elements which 
are not, as yet, available as practical alloys, shows certain important general trends. 
Firstly, the higher the melting point, the higher, generally, is the density. This 
tendency is illustrated by Fig. 1, which also shows the general trend in elastic modulus. 
Secondly, there is a close relation between melting point and the activation energy for 
those metallurgical processes dependent on atomic mobility, such as diffusion or 
high-temperature creep (Fig. 2). The exact significance of this second point would 
require a fairly detailed metallurgical discussion but it does mean, generally speaking, 
that for any given set of experimental conditions the rate of the relevant process, such 
as creep, is slower in those materials with a higher activation energy. Thus much is to 
be gained, it would seem, by moving to materials with higher melting points, the 
greatest penalties being those of weight, oxidation resistance and, usually, fabric- 
ability. It is not necessary, in the present context, to pursue these questions in detail ; 
they have, anyway, been discussed elsewhere“). The main point is that, from a survey 
of the materials situation in this light, two elements emerge as strong non-conformists: 
beryllium and the graphitic form of carbon. Both materials have a low density 
(beryllium 1-82 and graphite approximately 2-0); they both show limited ductility at 
temperatures low compared with their melting points—at room temperature, 
beryllium has a ductility of less than one per cent and graphite a zero ductility. 
Both have hexagonal-type structures, but with important differences, and both are 
materials of serious significance to advanced engineering. An important distinction 
between them at the present time is that, while beryllium is extremely expensive, 
graphite is relatively cheap—at least in conventional form. Even the newer types of 
graphite, which are certainly technically superior to those manufactured in the 
traditional manner, are unlikely to be ruled out on grounds of cost. There is also a 
very great distinction between the two in the matter of temperature resistance: 
indeed, graphite might be regarded, from one point of view, as the ultimate high- 
temperature material. The energy required to remove a carbon atom from a layer 
plane in the graphite lattice is very high—higher than that required for any other 
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Figure |. The trend of the room-temperature modulus and density of the elements with 

melting point. Beryllium and graphite are distinct exceptions to the general trend. Graphite 

(indicated by C*) has no melting point at atmospheric pressure, sublimating directly to the 
vapour phase at approximately 3,600-3,700°C. 


known lattice—and consequently, as graphite sublimes to the vapour phase without 
any intermediate liquid phase at normal pressures (at about 3700°C.), very large 
amounts of energy can be removed by ablation. This is of very great advantage, as 
those materials which melt are liable to be carried away from the surface in the liquid 
phase, and only the latent heat of melting is utilised. Unfortunately there are factors 
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Ficurr 2. Activation energies for creep and for self-diffusion plotted against melting point. 


which make it difficult (at least at the present time) to take full advantage of the 
possibilities that graphite theoretically offers. An assessment of these necessarily 
requires some appreciation of the fundamental nature of graphite. This is now 
considered. 


2. The Fundamental Nature of Graphite 


Graphite is a particular form of carbon structure, the atoms being arranged in 
planes made up of carbon rings, the carbon atoms in the rings being spaced at a 
distance of 1-45 A. In the more common form of graphite, these planes are stacked in 
the ABABA . . . . sequence, that is, every other plane is in alignment, as Fig. 3 shows. 
The spacing between the planes is 3-348 A and the interplanar bonding is fulfilled by 
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FiGuRE 3. The crystal structure of graphite. 


van der Waal’s forces. In the hexagonal metals, each layer plane is close-packed 
(each atom has six nearest neighbours, which is equivalent geometrically to adding 
an extra atom to the centre of each carbon ring) and these planes are arranged in the 
stacking sequence ABABA ... . as with the common hexagonal modification of 
graphite. There is X-ray evidence for the existence of a rhombohedral form of 
graphite, which is identical in its layer plane structure, the difference being that the 
planes are stacked in the ABCABC . . . . sequence. This is illustrated by Fig. 3. 
Such a stacking sequence is characteristic of the face-centred cubic metals. The 
rhombohedral form is unstable at high temperatures, and can be annealed, as it were, 
to the hexagonal form at temperatures above 1300°C., at least in artificial graphites. 
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Like many metals, graphite exhibits twinning, which can be initiated fairly 
easily. 


Graphite certainly differs in a number of significant ways from metals: the 
character of its interplanar bonding is an obvious example of such a difference. It is 
also distinguished by the fact that intermediate types of structure, between what 
might be called amorphous carbon and perfect graphite, exist and indeed represent 
stages in the graphitisation processes. Such structures are probably best described 
as ones in which a high degree of order within planes has been achieved—that is, 
carbon ring formation—without any highly developed ordering in their stacking. 
This type of arrangement has been termed “‘turbostratic”. These certainly represent 
important differences, because they do mean that conventional graphites may (and 
almost invariably do) contain regions which are highly disordered. Such regions have 
no counterpart in metals: the nearest parallel would be the structure of grain 
boundaries, but this is not a valid comparison, as such structures are merely highly 
localised atomic accommodations between two grain orientations. Indeed, such 
boundaries also exist in graphite. The fact is, then, that materials which are commonly 
described as graphites may really be no more than highly graphitic, containing regions 
which are not graphitic at all, with, in addition, considerable flaws and holes. 


Another important difference (with real practical applications) is that the graphite 
lattice is more anisotropic than that of any known metal. This anisotropy covers 
mechanical, thermal and electrical properties, and arises, of course, from the quite 
different types of bonding within and between the layers. 


In spite of these differences, however, the graphite structure, as such, is one which 
is amenable to consideration in terms of metal physics. Indeed, at the present time a 
very great deal can be gained by adapting dislocation theory to the case of the graphite 
lattice. Such a theory can go a long way, for example, in accounting for the twinning 
and stacking fault observations®: »), and developments in this direction will probably 
prove of high importance in elucidating the mechanical properties of graphite, as well 
as the kinetics of the graphitisation process itself. The situation might be summed up 
by saying that graphite should be accepted as a metal with polymeric tendencies. 


3. Graphite Technology 


Details of methods of graphite production need hardly be described in a paper of 
this type. However, some consideration must be given to the manufacturing questions 
if advances in graphite technology are to be given any just assessment. 


First of all, graphite does exist as a natural mineral and often relatively large and 
pure pieces can be mined—for example in Ceylon. Deposits exist in, among other 
places, Madagascar, Mexico, Korea, China and Germany, each having its own 
characteristics and each being utilised, in one way or another, by the graphite 
industry. As they stand, the natural products are of no commercial use, and must be 
converted into useful compacts. This can be done by mixing with a suitable binder 
material and graphitising the whole aggregate or, in one special case, actually com- 
pacting small graphitic particles, from which all impurities have been removed, by 
pressure at high temperatures. The resulting material is very highly oriented and has a 
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high density of about 2-1, compared with the theoretical density of the graphite 
lattice of 2-26. This figure is higher than that for the more common artificial 
graphites, which usually have a density of about 1-8 - 1-9. 

Natural graphite can be a useful addition to industrial mixes, but is not in itself 
the primary source of commercial material. Briefly, the essential production process 
is as follows. 


Carbon artefacts are prepared from the products of the distillation of petroleum 
and coal. The hard cokes derived from petroleum residues or from the pitch of the 
coal-distillation process are crushed to powder form, calcined to remove the volatile 
constituents, and graded into powders with suitable particle size distributions. Lamp 
black, carbon black or natural graphite fines may, or may not, be blended in. A 
binder material, which is usually ground pitch, is then added in suitable proportions 
and the whole is mixed. The mixture is then pressed or extruded. These “green” 
mixes are then baked at between 900 and 1300°C. under non-oxidising conditions, 
with the object of achieving the maximum carbon with the minimum of mechanical 
damage. The blocks of carbon so prepared are then packed into carbon powder and 
heated to not less than 2500°C. by direct current heating, during which the process of 
graphitisation proper occurs. 

This process of graphitisation is one which has been subjected to much research, 
but it cannot be said that any real understanding of the controlling factors has yet been 
achieved, even though a fair amount of knowledge of the character of ithe structural 
changes has been acquired™. As graphitisation proceeds, the linear dimensions of 
the graphite crystallites start at about 10-20 A, but increase rapidly for temperatures 
over about 1000°C. (which is the temperature at which a rapid evolution of hydrogen 
occurs). Growth, therefore, seems to be associated with the re-combination of edge 
atoms made available for bonding by the removal of their attached hydrogen or 
hydrocarbon groups. For carbon black, the dimensions parallel to the planes, L,, 
may reach 100 A or more (i.e. about 40 rings), depending on the particular material 
and that normal to the planes, L., may reach 50 A or more (i.e. about 15 planes). 
This alignment occurs at first without any rotational adjustment, so that orientations 
about the plane normals are random. As the temperature is raised, pairs of neigh- 
bouring planes assume the graphite relation, the actual temperature at which this 
occurs depending again on the specific material—it can be over 2800°C, although 
usually rotations begin at about 2000°C. There does seem to be some correlation with 
L, in this respect, as nearest neighbour graphitic ordering appears to start when L. 
reaches about 100 A. It has also been possible to detect second neighbour ordering 
and this has shown that, in the early stages, groups in the ABA and ABC sequences are 
of equal probability. The retention of bands in the ABCABC . . . . sequence would, 
of course, result in rhombohedral stacking faults. For prolonged heating, more and 
more layers orient themselves into the graphitic sequence. 


The graphite structure is obviously anisotropic and the coefficient of thermal 
expansion in the c-direction (normal to the layer planes) is about thirty times greater 
than in the a-direction. In cooling, after graphitisation, very great internal stresses can 
theoretically arise and these could be large enough to break the bonds and form voids. 
It is argued that the restraints imposed on polycrystalline graphite in this respect are 
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largely responsible for the inability to achieve the “perfect” graphite density of 2-26, 
although, of course, the effect may be mitigated, or obscured, by the fact that the 
starting compacts already contain voids which may not be filled by the graphitisation 
treatment. This could lead to anomalies. However, it may be important to recognise, 
in this connection, that the idea of void-creation by anisotropic shrinkage is tenable 
only if no stress-relieving (recovery) mechanisms operate. In particular, if dislocations 
are accepted, then the necessary geometrical adaptions during cooling could be readily 
provided by dislocation generation and movement. It seems likely that some such 
adaptation does, in fact, take place. 


The conventional methods of graphite manufacture necessarily involve the 
creation of micropores in the material. Considerable research has been devoted, 
in more recent times, to producing impermeable graphites, particularly for nuclear 
applications. It might be convenient to divide the pores into two broad types; those 
which arise at the particle boundaries in the baking processes, due to the loss of the 
volatile constituents, and those on a smaller scale which make up the misfit boundaries 
between adjoining crystallites of different relative orientation. While the former may, 
theoretically, be processed out, there is a limit to the latter contribution, set by the 
crystallite size. The larger the crystallites, and the greater the degree of preferred 
orientation, the less the microporosity. This has been confirmed by experiment, 
showing the crystallite boundary micropores to be about 3 to 8 A in width. This 
compares with interparticle voids of, say, 10,000 A or more. 


To sum up: the processes of graphitisation have certain general features which 
can be recognised and measured, but in detail they depend greatly on the materials 
used and on the exact graphitising conditions. This variability has made it difficult to 
establish a proper quantitative theory. The reasons why some carbons graphitise 
far more readily than others is not understood. For example, an anthracene low 
temperature char graphitises much more readily than chars from simple benzene 
derivatives, but, if the conditions are adjusted, graphitisable carbons can be prepared 
from benzene itself. Not only do the size of the crystallite parameters vary from one 
type of carbon to another, but so also do the spacings. 


It is possible to build up less and less defective graphite by polymerising smaller 
aromatic molecules at progressively higher temperatures. X-ray studies show that 
the progress of polymer formation takes place more rapidly with some than with 
others: carbon polymer from polyvinyl chloride gives a density of 2-25, whereas 
from polyvinylidene chloride the density limit is 1-59, even after 3000°C. heating. 


It does seem that certain types of crystallite boundary are extremely stable, even 
at very high temperatures, so that if the initial structure is such that these types are 
developed, then further growth ceases. This is clearly true for lamp black particles 
where the crystallites grow to the maximum size consistent with the retention of 
layer-plane surfaces. Surface planes which contain edge-atoms of the layer planes are 
therefore of high-energy. Growth which involves their creation is therefore impeded 
unless the energy can be derived from, say, an internal stress field. 


The processes of graphite manufacture described here have been evolved to meet 
the demand for a cheap material for use in electrical and chemical plant. There has 
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been no commercial incentive, until recent times, for any very high quality, and 
possibly expensive, product. The situation was changed, in the first place, by the need 
for an impermeable nuclear-grade material and this stimulated quite new researches, 
many of which were undertaken, and developed successfully, outside the established 
graphite industry. The most important of these may be classed as “impregnated 
graphites’’, which are considered separately in the following section. It may be 
observed here that a new stimulus to graphite research and development is now 
afforded by high-temperature technology, particularly, at present, by rocketry and by 
re-entry vehicles. Unfortunately, the market for such applications is negligible 
compared with the huge demands of industrial plant, so that the graphite industry has 
little incentive to undertake the expensive research and development work necessary 
if a radical change in the quality and properties of graphite is to be achieved. Much 
fundamental work is in hand within the industry and this is of the highest quality, 
but the results are not freely available and, certainly in matters in any way relevant to 
the processes of production, the closest secrecy is maintained. In an industry where 
success depends upon specialised techniques, developed to their present state by 
prolonged experience, such caution is understandable, but it may well be questioned 
whether it is, in the long term, altogether wise. For one of the consequences of this 
impenetrable barrier which the industry has set up around itself has been to compel 
groups with special interests to undertake, on their own account, new lines of work 
which have yielded advances of great consequence and promise. It can hardly be 
said to be a credit to the established graphite industry, with its unparalleled resources 
for research in the field, that the two major advances in graphite technology over the 
past few years have sprung into prominence because of work carried out in organis- 
ations external to that industry. These two particular developments concern impreg- 
nation techniques, and the consequent achievement of highly impermeable materials, 
and pyrolytic graphites, produced by the thermal breakdown of gaseous compounds, 
such as benzene, methane and propane. This situation is even less understandable 
when it is realised that the graphite industry has conducted successful research in these 
fields for many years and has certainly been in a position to launch out with new 
ventures, but, presumably for economic reasons, has not done so. The result is that at 
present the materials which seem most promising for aeronautical and space projects 
are being developed in newly-created organisations which lack the enormous back- 
ground of knowledge and experience which the major graphite companies have 
acquired over the years. 


4. Pyrolytic Graphite 


Carbon deposited from a suitable vapour at sufficiently high temperatures (say 
2000°C. and above) can be highly graphitic in structure, with negligible permeability, 
and strongly oriented crystallographically, the layer planes lying parallel to the 
deposition surface”). Such graphite sheets, termed “pyrolytic”, are therefore 
extremely anisotropic, with low conductivity normal to the plane of the sheet and a 
high conductivity in the plane. Some of the principal properties of products now 
commercially available are listed in Table I. 
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The successful development of techniques for producing large, shaped pieces of 
pyrolytic graphite has been one of the noteworthy advances of the past year or so. 
Control of all the process variables, to the required degree, is not easy. For successful 
products, the conditions of deposition should be constant throughout, which means 
that the surface at which the hydrocarbons are being broken down and deposited 
should be held at constant temperature. As the deposit grows, it acts as a layer of low 
thermal! conductivity between the wall and the deposition surface (usually being 
oriented so that the c-direction is normal to the growth face) and this upsets the 
thermal distribution. The cross sections of the deposits, when examined micro- 
scopically, reveal a cone-like growth from the base material. The nature of the 
deposit is, incidentally, affected by the nature of the substrate, ceramic bases being 
coated with a compact carbon, while metals and metal oxides give a powdery deposit. 
It does now seem that fine-cone structures are associated with a high resistance to what 
is called delamination, or splitting along the layer planes, and the deliberate attempts 
to produce such materials has resulted in a striking improvement. It has also been a 
commonplace assumption that higher graphite strength is necessarily associated with a 
higher density material. This does not seem to be well borne out by some of the 
pyrolytic work. It is possible to increase the density of such graphites by high tem- 
perature annealing, but the resultant strength can be less by a factor of two or so. 


5. Impregnation—Filling 

The porosity which exists in conventional fine-grain artificial graphites takes the 
form of interconnecting macropores, about 1 in radius, forming a network 
between the original particles. The possibility of improving the properties of standard 
commercial grade graphites by attempting to fill these pores is no new idea. It is not, 
however, so easy to achieve a radical improvement by this means as might be thought. 
The character of the initial graphite is important, and so also are the details of the 
impregnation process. Obviously the main interest in this development has been in 
connection with the use of impermeable graphite in nuclear engineering®?. The 
improvements in physical properties generally are of significance, however, in other 
applications and Fig. 4 shows the extent of these improvements compared with the 
original starting material. 


Technically, the long periods of treatment necessary for the production of 
impregnated graphites are a disadvantage, but if it were proved that these graphites 
were particularly outstanding for nozzle applications, then the time factor might well 
be accepted. It is understood that some new methods of manufacture are now being 
attempted by gas-cracking in beds of carbonised matted fibres. 


The principal method of liquid impregnation depends upon the use of furfuryl 
alcohol (C;H«O2) plus a suitable catalyst, so that polymerisation occurs, yielding a 
solid resin which is carbonised at 1000°C. to give about 50 per cent carbon. During the 
carbonisation the resin shrinks uniformly without becoming liquid. As the imper- 
meability is unaffected by high-temperature treatment, the method may be applied to 
carbons which can subsequently be graphitised. Fine-grained materials appear most 
suited to this kind of treatment, and it is usual to carry out two impregnations. In 
the case of tubes, a common one for reactor applications, one impregnation is made 
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Ficure 4. The degree of improvement in the properties of conventional graphites brought 
about by impregnation. 


from each side. Considerable improvements in mechanical properties, such as strength 
and modulus, result from this kind of treatment, as indicated by Fig. 4. The successful 
use of furfuryl alcohol has naturally led to its being considered as a binder material. 
Quite large section bars of graphite have been produced with such a binder, and the 
permeability of the product normal to the extrusion direction was fifty times lower 
than that of a conventional low-permeability graphite. 


Impregnation can also be carried out using normal pitch binder material, but a 
large number of treatments is required to change the permeability by a factor of a 
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thousand or so and the process of decomposition is in any case necessarily slow. 
Impregnation by a concentrated solution (7 per cent by weight) of sugar has proved 
highly effective. On heating slowly to about 1000°C., the sugar decomposes to leave a 
residue of carbon in the pores; four to five impregnations can reduce the per- 
meability by a factor of about 106. About six or seven such impregnations, followed 
by one or two impregnations of furfuryl alcohol, result in an approximate doubling of 
both modulus and strength. 


It is also possible to carry out what is virtually a gas impregnation technique, 
although this might well be classified as a pyrolytic deposition method, particularly as 
coatings can be applied successfully by this means. In the particular impregnation 
technique of most interest the pyrolysis is so controlled that it occurs only within the 
pores of the material, the sealing-up action therefore takes place wholly internally 
and no surface coating problems arise. With a nitrogen-benzene vapour atmosphere, 
the transition from surface to internal deposition occurs between 750°C. and 800°C. 
The weight can increase by about 9 per cent without any dimensional changes, and 
permeabilities of 1-5 x 10-9 cm.2/sec. can be achieved. 


This temperature range applies to graphites with a mean pore size of about 0-5 
micron. A wider temperature range would be needed to make use of a coarser 
graphite of larger pore size. 

One factor of some interest to the high temperature field is that the oxidation rate 
of the finally graphitised products of this type is about one-twentieth of the original 
graphite, although these figures relate only to performance at a relatively low 
temperature (550°C.). 


6. Coatings 


Graphite is readily oxidised at high temperatures and for successful operation 
some form of protective coating is required, such as an impervious layer of 
carbide'®'!2), A large number of research organisations are now engaged in examining 
methods of coating graphite and, because of the commercial implications, few details 
can be ascertained. One method is based upon the vaporisation of halides (NbCIs, 
Ta Br;, Zr Brs, for example) at temperatures in the region of 2000°C., at low pressure— 
a few centimetres of mercury. The carbide coatings deposited (NbC, TaC, and ZrC in 
the cases quoted) bond well to the graphite, being locked into the porous structure, 
and they can be made in thicknesses of the order of 1 mm. They appear resistant to 
thermal cycling and, by using the base carbon to form the carbide, the rate of growth 
is controlled by diffusion and is, therefore, largely insensitive to the vapour flow 
conditions. 


Profound disagreement exists at present regarding the improvements that the 
available coatings offer. In the case of silicon carbide (SiC), some quite striking 
improvements in resistance have been reported for 30 seconds exposures to a stag- 
nation temperature of 3000°K: the actual surface temperature of the carbide may not 
rise above 2000°C. Other workers in the field have expressed the view that this is an 
optimistic estiinate and, if anything, a coating of 60W 25Cr 15A1,0; is superior in this 
temperature range. At the lower temperatures, improvements are easier to obtain 
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using coatings such as silicon nitride (SiN) and boron nitride (BN). Boron nitride 
oxidises at 800°C. and silicon nitride oxidises rapidly above 1400°C. Expert views in 
the matter of coatings research differ very sharply. Many believe that the main weight 
of research should be concentrated on developing graphite itself to its most resistant 
form. There is now some sign that this view might be well founded in that pyrolytic 
graphite has advanced technically very markedly over the past year. Another view 
is that admixtures of carbides with graphite will prove far more profitable than 
coatings, as such. Finally, another body of opinion feels that the whole science of 
coatings is in its very early stages and that continued research in this direction will 
ultimately pay high dividends. 

There seems no reason at present to assume that improvements along each of 
these three lines could not be exploited to advantage by future technology. There is 
little doubt that the oxidation resistance of graphite is an important factor, even 
though for many applications the atmosphere is reducing and not oxidising. If 
research were to be seriously influenced by this last consideration, then the results 
might be unfortunate, as it may well be necessary in the future to take up the subject 
again if highly oxidising fuels are brought into use. 


7. Graphite as a Future Structural Material 

The foregoing survey has drawn attention to some of the special features which 
distinguish graphite from other materials, and has indicated the principal classes of 
graphite at present available: these may be grouped as (a) natural, (b) conventional 
artificial, (c) impregnated and (d) pyrolytic. 

If any just assessment is to be made of the part that graphite might play in future 
technology, then some special examination is required of its mechanical 
properties *!9), 

The ductility of graphite at temperatures up to about 1500°C. may be taken to be 
zero. At 2500°C. it is about 5 per cent; at 2700°C., 10 per cent, and about 20 per cent 
at 3000°C. The basic mechanics of the crystal lattice of graphite, and the exact role 
that boundaries play in deformation, are so little understood at present that it would 
be unwise to draw the conclusion that ductility is an inescapable limitation. The stage 
has now been reached in graphite studies where a real advance in understanding is 
possible, rather on the lines of that which has taken place in metal physics over the 
past twenty years or so. It should be remembered, however, that the progress of 
fundamental studies in crystal mechanics has been made possible, to a large degree, 
by measurements and observations on the behaviour of single, isolated, metal crystals, 
and on limited aggregates of crystals. Unfortunately, the successful production of 
graphite crystals, large enough and perfect enough for serious experiment, has not yet 
been achieved. The importance of this course of fundamental research cannot be 
over-emphasised. 

One of the interesting facts about the conventional graphites is that their strengths 
and elastic moduli increase with temperature—at least up to a limit, above which they 
fall off (Fig. 5). The tensile strength, for example, can approximately double from, 
say, 3 x 105 Ib./in.2 at room temperature to 6 « 103 Ib./in.2 at about 2500°C. The reason 
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Figure 5. The variation of the strength of graphite with temperature compared with that 
of other materials. 


for this increase is not established. It is tempting to think that it is a secondary effect, 
due to the enhanced plasticity at the higher temperatures, which in turn reduces the 
severity of any stress concentrations present. In other words, the ability of the 
material to work-harden at the higher temperatures more than compensates for any 
true thermal weakening. The same sort of effect is, indeed, observed with other 
substances, such as the ceramic oxides. It might be expected, then, that the more 
perfect the graphite, the less evident this strength increase with temperature should be. 
Such an expectation is not borne out by experiments on pyrolytic graphites, which are 
highly oriented, and more perfect, and show an increase in strength from about 
18 x 103 Ib./in.2 at room temperature measured parallel to the highly oriented layer 
planes (this figure being six times higher than that for conventional graphite) to 
about 33 x 103 Ib./in.2 at 2500°C. and 50 105 Ib./in.2 at 2700°C. This strength at 
2700°C. puts graphite into a class of its own, far higher than the strengths of the 
refractory metals, rhenium, molybdenum, tantalum and niobium. Few tests have, 
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in fact, been reported on metals at this temperature, but even at 2000°C. the highest 
metal strengths, under short-duration loading, are less than 10 105 Ib./in.2 and 
there is no recorded strength as high, even, as 5 x 103 Ib./in.2 at 2700°C. The strength 
of pyrolytic graphite is, then, at least an order of magnitude better at these extreme 
temperatures. It should be remembered that such graphites are highly anisotropic, 
and the strength perpendicular to the layer planes is much lower. Although, then, 
more perfect types of graphite still show the temperature-strengthening effect, the 
original idea that this results in some way from increased plasticity need not be dis- 
carded outright. It does mean that gross inhomogeneities, macropores and the like, 
are not responsible, and it must be supposed that there are features in the lattice itself 
which can confer brittleness at low temperatures while being decreasingly effective 
in this respect as the temperature is raised. The problem is therefore likely to involve 
dislocations and much must await the outcome of studies in this direction. 


The increase in the elastic modulus of graphite with temperature is also a useful 
feature as far as design is concerned, but is even less well understood than the strength 
increase. Again, the amount of the change in modulus is a function of the direction in 
which it is measured—that is, parallel or perpendicular to the layer planes in pyrolytic 
graphite, or, with conventional graphite, parallel or perpendicular to the direction of 
extrusion. With pyrolytic graphite (“‘parallel’’ measurements) the modulus increases 
from about | -0 x 10° Ib./in.2 to 1-5 « 10¢ Ib./in.? at 600°C. 


Conventional graphites, which have about the same Young’s modulus parallel 
to the extrusion direction, increase by about the same amount over the range 0-2000°C. 
The rigidity modulus of such graphites also increases by nearly 50 per cent up to about 
1500-2000°C., above which it falls off. 

The creep of graphite exhibits much the same general features as that of metals 
and appears to be capable of a physical treatment on similar lines®®. Some creep 
curves are illustrated by Fig. 6 and an examination of this type of data in terms of 
activation energy (that is, by plotting creep rate on a logarithmic scale against the 
reciprocal of absolute temperature for a constant stress) leads to results of the kind 
shown in Fig. 7. Without entering into any detailed discussion of the physical factors 
which affect the creep behaviour, it is perhaps worth noting that the values of creep 
activation energy can be very high—around 200 kcal./mole, approximating to the 
energy for self-diffusion. Graphite is, then, relatively creep-resistant, in the sense 
that a high energy is required to activate the elementary creep process. Some caution 
is required in interpreting this factor, for a high activation energy alone does not 
necessarily imply a low absolute creep rate. However, in general, resistance to creep 
at any given temperature is strongly related to the activation energy, and high values 
are certainly desirable. 

Another factor related to that of the creep activation energy is the energy 
associated with ablation. As was noted previously, the sublimation energy (for the 
direct solid-vapour transition) of graphite is high, and consequently absorbs a large 
amount of any heat developed aerodynamically, provided that this heating is 
sufficient to promote the process in the first place. Depending on the material, 
actual burning, or depolymerisation, or melting, are all, of course, processes capable 
of contributing to this energy absorption, but none is theoretically as efficient as the 
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Ficure 6. Creep tests on graphite at 2,000 lb./in.? parallel and perpendicular to the grain 
at the following temperatures : — 
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sublimation of graphite, which is capable of absorbing about 20,000 B.t.u./Ib., 
compared with figures of about 5,000 B.t.u./lb. for plastics, 1,000 B.t.u./Ib. for 
beryllium and 200 B.t.u./Ib. for copper. Unfortunately, conventional graphite is a 
relatively good conductor, with an average coefficient of thermal conductivity, at the 
lower temperatures, equivalent to magnesium or beryllium. This decreases very 
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Ficure 7. Creep rate plotted against 1/7, as derived from the curves of Fig. 6. 


usefully as the temperature is raised, so that at 800°C. the conductivity in the low- 
conductivity direction (perpendicular to the layer planes) is down to about one-tenth 
of its value at room temperature. At temperatures of 2000-3000°C. the value is 
probably as low as one-hundredth of the room-temperature conductivity. In pyrolytic 
graphite, formed as, say, a casing material, with the layer planes parallel to the 
surface, this low conductivity will usefully limit the transmission of heat to the 
interior. The conductivity parallel to the surface for such graphites is higher by a 
factor of over 100. This anisotropy in the conductivity may not be of critical import- 
ance in thermal protection systems, but raises special design problems for more 
general structures. The anisotropy in thermal expansion is more troublesome. In the 
approximate temperature range 0-200°C., the thermal expansion coefficients differ 
by a factor of 100 in the “parallel” and “perpendicular” directions; at higher tem- 
peratures, the ratio approaches 10. Expansion is greater perpendicular to the layer 
planes (25-0 10-9/°C.) than it is in the parallel direction (~0-2 10~-6/°C. at 
150°C., ~2-0x 10-¢/°C. at 1000°C.), and if the graphite were deposited on a metal 
substrate, the differential expansion parallel to the surface could be serious, steel and 
beryllium being about 12x10-¢/°C., and aluminium and magnesium around 
25x 10-9/°C. Tungsten, however, is usefully low—about 4x 10-9/°C. It should 
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also be noted that similar differential expansion problems would arise if pyrolytic 
graphite were coated on a conventional graphite base. Indeed, experiments on such 
materials in plasma jets have shown that the surface deposit can be fractured by the 
thermal stresses so created. If, however, the degree of orientation in the graphite 
could be varied with depth, the surface layers being highly oriented (layer planes 
parallel to the surface) giving a low conductivity in the perpendicular direction, and 
with the material becoming gradually more isotropic in the interior, then real benefits 
might be achieved. 


Although high-purity graphite is an abrasive as far as cutting tools are con- 
cerned, it can be machined fairly easily, and this is the usual way of manufacturing 
components. Graphite-graphite bonding is possible, using either a conventional 
binder or a substance such as furfuryl alcohol, and finally re-graphitising, or by 
brazing techniques®”. Copper-silver eutectic alloys have been used in the latter 
method, by heating in vacuo at 950°C., as have titanium and zirconium heated in 
vacuo to 1850°C. Metals can also be joined to graphite by some such brazing tech- 
nique, so that complex, bonded assemblies are certainly possible, even if technical 
difficulties still remain. 


With such evident virtues it might well be asked why it is that graphite has not 
been very seriously considered as a high-temperature structural material, even though 
it is now accepted as the best material for certain specific applications—for example, 
rocket nozzles. The main reasons are, firstly, that its properties at low temperatures 
are not acceptable for primary structures, in particular its negligible ductility. The 
pyrolytic graphites are no better in this respect, although they are somewhat superior”? 
above 1650°C. Secondly, conventional graphites show a notorious variation in their 
properties and no kind of specification exists at the present time. This situation will 
certainly be improved by the development of the newer types. Thirdly, and perhaps 
most significant of all. graphite reacts with gases at high temperature and, in particular, 
burns away in an oxygen environment. If the gas stream impinges on the surface at a 
high velocity, then erosion also occurs and it is this combined process which really 
causes most concern at present. Very little is really understood about the mechanics 
of erosion itself, so that the part that the particular characteristics of the chemical 
attack might play in the case of graphite are even more difficult to assess. Much 
fundamental work is now going on, and it is too early to attempt any prediction of the 
way in which a radical practical improvement might be brought about, but, as this is 
certainly a major question, a special review will be made of the underlying factors in 
the following section. Before doing this, it may be timely to comment on the 
frequently-expressed view that, as the environment in rocket nozzles is usually 
reducing, not oxidising, the attempts to decrease the oxidation rate at the surface are 
mis-guided. The main purpose of fundamental surface studies is to elucidate the 
basic mechanisms of the reactions, and oxidation, apart from its intrinsic interest, 
provides a clue to some of the more general characteristics of such reactions in 
graphite. In any case, it would surely be short-sighted to abandon the search for a 
means of reducing oxidation merely because, at this particular time, in the evolution 
of a particular technology, oxidation is not the primary problem. 
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8. Surface Reactions: Some Fundamental and Practical Factors 


In the conditions of the rocket nozzle, graphite is eroded by the burning gases. 
This is particularly severe with the newer high-temperature fuels, such as those 
containing aluminium. Opinion on the exact mechanism of this erosion, and the 
factors influencing it, is divided. A microscopic study of certain graphites, which have 
been widely used in nozzles, after exposure to rocket conditions, has resulted in the 
conclusion that the attack occurs preferentially at the interfaces of the original 
particles—that is, in the binder material. Another view is that the presence of metallic 
particles may seriously accelerate the wear by abrasive action. Yet another is that the 
effect is worse with the metallic propellants, simply because they operate at higher 
temperatures and for no other reason. This view is almost equally counterbalanced 
by the opposite, evidence being cited for the high resistance to erosion of graphite 
under the action of high-aluminium fuels. This dissension really underlines the fact 
that erosion is a chemically-activated process and not a question of wear mechanics, 
at least in the early stages, and that it is very sensitive to the exact combustion 
conditions. 


The possibilities of silicon carbide-silicon boride coatings have been previously 
noted, attention being drawn to the disagreements which exist at present regarding the 
temperature range over which such coatings may prove beneficial. In fact the non- 
uniformity of the attack, once the carbide coating has been eroded through locally, 
results in a distortion of the nozzle profile which is more deleterious than that brought 
about by a more uniform wear. Silicon carbide coated graphites fail by the con- 
version of the silicon carbide to silicon dioxide, which liquifies at 1600°C., and is 
driven off by the gas stream. 


As with all chemical kinetics, the rate of attack, other things being equal, will 
increase with temperature. The temperature factor alone, however, cannot explain 
the wear rate, and this explanation may be discarded at the outset. The constituents 
of the fuels are, however, of great importance, as has been demonstrated by experi- 
ments in which magnesium and aluminium fuels were used in similar nozzles at very 
nearly the same temperature. Only the aluminium-containing fuel eroded the graphite 
to any serious extent, and this suggests that the important factor is a chemical reaction 
between some constituent, such as Al20;, and the graphite. 


Fundamental studies are clearly required for a proper assessment and certain 
results recently presented on work which is still in progress °) are drawn upon here. 
The nature of the chemical attack on graphite surfaces in reactions with CO2, O2, Os, 
and O, at temperatures below 1000°C. has been studied by optical and electron 
microscopy. Hennig finds that metallic particles of Fe, Au, Ag, and W do not act 
individually in catalysing the reactions, but only when aggregated into larger colloidal 
particles of about 200 A diameter. The attack is not on the basal plane but only on 
the edges of the steps, and then only in certain preferred directions, so that a sort of 
channelling results. If vacancies are present, however (as produced by irradiation, 
for example), then the laver planes are attacked at once and pitting results—presumably 
this accords with the previous observation, the attack being at the edges of the 
vacancy clusters. Whatever the detailed interpretation of these and other related 
effects may be, it is clear that the crystallographic nature of the surface, and the 
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constituents of the gases, are of fundamental importance. The most beneficial surface 
structure appears to be that presented by the layer plane face, at least for the case of 
O>; in fact, chemisorbed O>2 protects the layer planes. On the other hand, O; attacks 
the layer planes direct. The presence of chlorine, bromine and carbon monoxide 
affects the nature of the attack and slows it down. Grinding, on the other hand, 
accelerates chemisorption, whereas subsequent annealing retards it. 


It is very likely that the binder material (or, more properly, the boundaries) offers 
sites for preferred attack by O2, because such boundaries may have a locally more 
disoriented structure, with a higher concentration of faults and vacancies, and with 
exposed layer edges. An aspect of this dependence has been demonstrated by natural 
graphite, of micron size, moulded without a binder up to densities of 2-0. Heat 
treatment at temperatures up to 3000°C. results in a gradual decrease in the re- 
activity, at a given temperature, to carbon dioxide; the factor is about 30. The 
activation energy remains the same (42 kcal./mole), so there is no question of a 
differing mechanism, the result being due to the reduction in the extent of the vulner- 
able material by processes of growth and boundary diffusion, and with the consequent 
dissolution of stress concentrations built up in the moulding process. 


The chemical considerations alone are still inadequate, however, in judging the 
behaviour when the erosion is superimposed upon conditions of thermal shock. 
Differential expansions in the surface may readily result in the development of micro- 
features which are susceptible to attack: the creation of surface steps by slip is an 
example. This is a field where the crystal mechanics of graphite are relevant, and 


this part of the fundamental theory is still in its infancy. The actual surface chemistry 
of graphite is indeed much better developed and paramagnetic resonance techniques 
have enabled a number of theoretical possibilities to be examined. The details of this 
work would lead into special technicalities which cannot be pursued here but it should 
be recorded that these researches are of high importance, and of great relevance to the 
high-temperature field. As with all properties, the rate of attack and its nature varies 
greatly between graphites of different origin and pre-treatment. This is consistent 
with the view that such attack is strongly related to the microstructural detail. 


Graphite has a good resistance to thermal shock, a coarse-grained graphite being 
better in this respect than the finer grades, which usually have better strength 
characteristics. Decreasing the amount of filler component decreases the thermal 
shock resistance, so that it would appear that the boundary material is largely 
responsible for absorbing the anisotropic thermal stresses. 


The degree of oxidation-pitting that occurs is dependent on the local con- 
centration of impurity in the graphite, being worse for titanium, vanadium and 
calcium than with iron and other impurities, at least up to about 1000°C. Va (giving 
Va20s) plus NaF is a particularly serious catalyst, and the attack (say 2-5 p.p.m. of 
each) is worse than the sum of the individual rates of each separately. The purity 
of the graphite is therefore important. 


The rate of reaction of carbon with gases is influenced also by the porosity and the 
amount of internal surface material, and by the diffusivity of the gaseous reactants. 
At low temperatures the dominant reaction is simply that of chemical attack at all 
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surfaces, internal and external, and the slower the rate the lower the concentration 
gradients. This would apply to the phenomena already discussed. As the temperature 
increases, the reactivity rate increases, the diffusion of the gas through the pores 
becomes a controlling factor, the effective rate decreases, and the overall activation 
energy, as measured, decreases. For even higher temperatures, the mass transfer 
through the boundary gas layer becomes important, this layer decreasing in thickness 
with increasing gas velocity. It is therefore to be expected that porosity is of impor- 
tance as a rate-determining factor in oxidation, but this is not significant for tempera- 
tures below about 600°C. or so. 


Returning to the nature of the burning of graphite, it has been observed that very 
special surface features are developed. All the layers of the steps which are attacked 
do not burn off at the same rate and, as a result, a sort of undercutting occurs because 
of the fact that some planes burn more slowly—these are spaced at about 0-5, with a 
very regular periodicity even in the most perfect graphite crystals obtainable: again, 
these experiments are at temperatures below 1000°C. 


In measurements made at very high temperatures (2400°C) with gas velocities 
in the region 20-60 metres/sec. in gases containing about 50 per cent H2O, there is 
some evidence that the rate is proportional to the square root of the velocity. The 
latest series of measurements indicate that the rate of reaction is in the region 3-5 to 
5-0 x 10-3 gm./cm.? sec. at gas velocities of about 50 metres/sec., this rate decreasing 
with graphite density. The difference between graphitised carbons and simple carbons 
appears to be small, at least at temperatures of about 2000°C. On the other hand, 
there is experimental evidence available which shows that pyrolytically deposited 
carbons have a significantly lower rate of reactivity than either commercial carbons 
or graphite, by a factor of about 3. In combustion tests of this kind small incandescent 
grains have been observed in the combustion gases and these appear to be eroded by 
the burn-up of the interparticle bonds. Incidentally, it is curious that the degree of 
parallel stacking appears to increase with oxidation; the effect being partly reversed 
when the oxygen is removed. No very convincing explanation of the effect has been 
deduced. 


There is an almost total lack of correlation between reactivity and many of the 
bulk measurements such as, resistivity, ash content, degree of crystallinity, so that 
only direct measurements have any real value. 


Pyrolytically-coated graphites have shown about the best resistance to erosion so 
far obtained but, at the very extreme temperatures, about 3000°C., they simply flake 
off if there is not a proper thermal matching to the under surface. 


Impregnation can be beneficial in a quite different way, in offering a means of 
energy extraction by vaporisation: for example impregnation with ammonium 
iodide can hold down the surface temperature for several minutes at a quite high rate 
of heating. Materials with higher boiling points are, however, required. Impreg- 
nation with metals appears unlikely to be of great value. Firstly, the rapid differential 
thermal expansion at the start of the combustion could result in the fragmentation of 
the surface layers and, secondly, the metal may well be swept out in the liquid phase, 
giving little benefit where it is most needed. For a reasonable quantity of such 
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coolant to be introduced, a fair pore size is required, so that the metal can be 
absorbed in the graphite structure—at least in the boundaries. 


Rare earth additions have been experimented with, such as cerium oxide and 
aluminium fluoride, but for operation above 2000°C. no really noteworthy advance 
has been achieved. One difficulty in drawing conclusions from the test results is that 
they differ greatly with the exact test conditions, so that materials which appear to 
show improvement in one jet torch are disappointing (and even worse than unimpreg- 
nated graphite) when tested in another. The possibilities of impregnation are, how- 
ever, scantily investigated, as yet, and some real advantages may yet be attained by 
such means. 


The actual admixture of carbides to the graphite mixes has also been investigated, 
with promising results. Combinations of graphite with BsC and SiC have resisted 
temperatures around 1600°C. for some 80 min. without failure by burning. No 
improvement has been achieved to date by metallic impregnation—for example, 
with copper. 

The actual sublimation of carbon itself, in the graphite form, is the most obvious 
mechanism to use to advantage. The difficulty is to arrange that such erosion is 
truly atomic, and not in the form of particles, as it most certainly is at present. 


9. Conclusions 

It will be evident from the foregoing that graphite offers certain quite unique 
possibilities as a material in high-temperature structures. These possiblities will not 
be fully realised unless a number of quite fundamental issues are better understood 


and unless those designing in graphite appreciate to the full its very special 
characteristics. The recent improvements achieved in graphites should encourage a 
certain optimism regarding its ultimate utilisation as a serious structural material, 
although in the immediate future there are several iuore limited applications where 
it could be of real value—as a turbine blade material, for example. In the longer 
term, it will certainly play an important part in the composite assemblies required to 
meet the specifications of advanced projects in the future. 
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Boundary Layers over Infinite Yawed Wings 


J. C. COOKE 


(Royal Aircraft Establishment, Farnborough) 


SumMaRY: A method of calculating turbulent boundary layers on_ infinite 
yawed wings is given, making use of a method of calculating turbulent bound- 
ary layers due to Spence and of an analogy between three-dimensional and 
axi-symmetric boundary layers. It is also shown that the displacement thickness 
is equal to that computed using chordwise components and that the streamwise 
momentum thickness is approximately equal to the chordwise momentum 
thickness. Shock-free flow and small boundary layer cross-flow are assumed. 


Introduction 


In dealing with the boundary layers of infinite yawed wings’: ”’ use has generally 
been made of the “ independence principle,” by which the chordwise flow may be 
calculated as though no spanwise flow were present. This only applies for 
incompressible or isothermal compressible laminar boundary layers. This paper 


gives a method of calculating the flow in the turbulent incompressible or com- 
pressible cases where the independence principle does not apply. It is based on an 
axially symmetric analogy” for general three-dimensional flow. 


It is first shown that, in spite of the fact that the independence principle does 
not always apply, the “ displacement thickness ” calculated by means of chordwise 
components is the true displacement thickness in the case of infinite yawed wings. 


Assuming small cross-flow in the boundary layer, which in effect means small 
incidence, the flow is calculated by the axially symmetric analogy, following a 
streamline of the external flow. For infinite yawed wings it is possible to modify the 
details of the analysis so that, in effect, the solution proceeds across the span normal 
to the leading edge. An equation for the momentum thickness is thus obtained and 
this is very little more complicated than the corresponding two-dimensional equation 
of Spence’. Other properties such as displacement thickness and streamwise skin 
friction follow. 


No attempt is made here to calculate the cross-flow component of flow. This is 
of little interest when one is assessing the effects of boundary layer on the properties 
of a wing whose incidence is low enough for the cross-flow to be small. This cross- 
flow component would of course be required to determine separation, but the method 
probably becomes less valid as separation is approached. 
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NOTATION 


a velocity of sound 

B defined in equation (12) 

c chord 

C, pressure coefficient 

c, skin friction coefficient 
F(M),G(M) defined in equation (12) 

h,,h, defined in Section 3 
H form parameter, =5,/6 


K=zT,./T.. 
M Mach number 
p pressure 


R_ Reynolds number 
r in Section 6, recovery factor =0-89 
r radius of analogous axi-symmetric body 
S$ arc measured along a streamline 
T temperature 
u,v velocity components along axes of x and y 
u’,v’ velocity components along axes of x’ and y’ 
v.=U, sin A 
U. velocity outside the boundary layer 
u.=U.? —U,,? sin? A 


x,¥,Z co-ordinates: x normal to leading edge along the surface, 
y¥ spanwise, Zz normal to surface 


, 


x,y,z’ co-ordinates: x’ along the external streamline, z’ normal to 


surface, Y normal to x’ and 2’ 
8B angle between x and x co-ordinate axes 
6, given by equation (1) 
8, displacement thickness given by equation (la) 
true displacement thickness 
6,6. streamwise momentum thickness 

=6(U.8/¥)"? 

A angle of sweep 

vy kinematic viscosity 

p density 

p defined in Section 4 
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skin friction 
¢(H),v(H),x(H) defined in Section 5 


Suffixes 
main stream values 


values at a reference station 
values at infinity 
recovery values 
values at a “mean” position 
incompressible values 
s values at an isentropic stagnation point 
T values at transition. 


2. The Independence Principle 


In laminar incompressible flow the independence principle is well known for 
infinite yawed wings. It is that the chordwise flow may be calculated as if no span- 
wise flow were present. It is sometimes thought that the independence principle may 
be stated as “ everything is independent of y,” where y is measured along the span, 
i.e. along the generators. It is important to realise that these are not the same. It is 
well known that the independence principle does not hold in turbulent or com- 
pressible flow’: » in general, yet if the wing were truly infinite one might expect the 
second statement “everything is independent of y” to hold even if the first does not. 
If steady flow is possible at all then the second principle must hold, and there seems 
to be no reason why steady flow should not be possible. It is assumed that the 
second statement, called here “ Principle II”, holds. 


3. Displacement Thickness 


In some circumstances it is necessary to know the displacement thickness 6* in 
order to assess the effect of the boundary layer on the external flow. 4* is defined 
to be the amount by which the body must be deformed or “thickened” in order that 
the potential flow about such a deformed body will give the same flow outside the 
boundary layer as the actual flow. We define 4, by 


Pelle 


oO 


u and u, being the chordwise components of the velocity inside and outside the 
boundary layer, p and p, the corresponding densities, and z the distance measured 
normal to the surface. 


We now show that the displacement thickness 6* is equal to 6,. The 
proof follows that of Lighthill for the general three-dimensional case. Lighthill® 
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Figure 1. FiGure 2. 


took an element of fluid in the streamwise direction and showed that, if x’ and y’ are 
measured along and perpendicular to the external streamlines (see Fig. 1), 


is not the true displacement thickness, because in obtaining it by flow through a 
control surface no allowance was made for the fluid flowing into the surface through 
its sides. When this was accounted for he found that the displacement thickness 
became 


pUhy ey pUdyhedx, 
Ud, oy’: 


where fA, and hd, are determined from the line element in the surface, that is 
ds? =h,?dx’? + h,*dy”. 


If axes x and y are taken normal to and along the generators, respectively, then 
the displacement thickness, as before, can be seen to be 4, together with a correction 
term analogous to that in equation (2). However, so long as everything is indepen- 
dent of y the correction term will vanish, since the fluid coming in on one side of the 
control surface is exactly balanced by the fluid going out the other side. An 
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analytical proof, starting from equation (2) and transforming the co-ordinates to 
chordwise and spanwise, can also be given. 


Hence, so long as Principle II holds, the displacement thickness 4* is equal to 
6,, aS defined in equation (1). This is so even in turbulent flow and compressible 
flow, for which the independence principle does not apply. 


4. Infinite Yawed Cylinders in Incompressible Turbulent Flow 


To investigate the development of the boundary layer over infinite yawed wings 
when the flow is turbulent, it is first necessary to calculate the laminar boundary 
layer up to transition. For this part, comprehensive accounts have been given by 
Rott and Crabtree“ and by Crabtree’. The Appendix gives an alternative method. 


After transition the independence principle no longer holds. Cooke'*’, however, 
has given reasons for believing that the flow in the streamwise direction behaves 
much as in two dimensions and, on that basis, has derived a momentum equation 
which applies if taken along an external streamline. He showed that, if 4 is the 
momentum thickness in the streamwise direction, and 


then, provided that the cross-flow is small, 


(eu. 
Os 


} -0-0106 


where U, is the external velocity along the streamlines, ds is an element of streamline 
length and p is a function determined by the external flow conditions. 


For an infinite yawed wing, Zaat®’ has shown that 


A? 


where A is a constant and 
=U,’ —U,,? sin’ A, 


U.,, being the velocity at infinity and .\ the angle of sweep. 
Hence as 


along a streamline. 


As long as Principle II holds, the integral along the curve A’ B’ (Fig. 2) can be 
replaced by the integral along AB, since any quantity which has a certain value at 
P’ has the same value at P. 
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In Fig. 2, dx=ds cos 8, where tan 8=7,/u., ve =U,, sin A, U2 =u? + 


U, Ox’ 
Integrating (3) with respect to x, and expressing U, u., 6 and x non-dimensionally 
in terms of U,,, and c, we have 


(4) 


Hence 


26"? = dx + constant, (5) 
where R=U.c/¥. 


In this equation x is distance along the surface of the wing normal to the 
generators. 


The equation depends on the assumption that the streamwise skin friction takes 
the same form as that used in two dimensions, namely 


=o.0088 

Doubts have sometimes been expressed as to the validity of applying this 
assumption to three-dimensional flow and these may possibly be justified in general. 
However the method has been applied“ to some unpublished experiments by 


DISPLACEMENT THICKNESS 
CALCULATED 
MEASURED 


0-002 
MOMENTUM THICKNESS 
CALCULATED 
© MEASURED 
0-001 
fe) | j 
0-4 0-5 O6 C7 0:8 o9 1-0 


Figure 3. Momentum and displacement thickness for swept wing R.A.E. 101, thickness /chord 
ratio 0-045, angle of sweep 55° (Brebner). 
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Brebner on the turbulent boundary layer at mid-semi-span of an untapered wing 
with a sweep angle of 55°. Fig. 3 shows reasonably good agreement between 
experiments and calculation by the present method for this wing and lends support 
for the assumptions made in deriving equation (5). 


Equation (5) will determine the streamwise momentum thickness at any point, 
starting from transition, assuming that its value there is known from the laminar 
calculations and that it is continuous at transition. 


Equation (5) may also be obtained by the use of the axisymmetric analogy as 
applied in Section 6. 


5. The Determination of H 


If it is required to find the displacement thickness in turbulent flow it is 
necessary to find H, where 


H=3,/6. 


Cooke“ found that, in an example resembling an infinite swept wing, H could 
be determined from the equation 


dH dU 


d. 
the same equation as used by Spence“ in two dimensions. In this equaiion 
(H)=0-00307 (H — 1)’, 
(H)=9-°525 (H — 1:21) 1). 


Integration of equation (6) gives 


5 


2 
U.2x (H) = constant — 0-00135 ds 


(-) 


0-442 
where x (H) =2-105 H-1 


Use of equation (4), as before, gives 


x 


3 
(H) =constant — 0:00135 dx. . (8) 


e 


The initial value of H, at transition, is usually taken as 1-4. If the flow is such 
that separation is nowhere approached it may be sufficiently accurate to take H as 
constantly equal to 1-5, instead of attempting to calculate its variation. 
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6. Infinite Yawed Cylinders in Compressible Turbulent Flow 


In the compressible flow case shock-free flow is assumed, except possibly at the 
leading edge. It is also assumed that » varies as the absolute temperature raised to 
the power 8/9. 


Spence’ has given a method of solution in two-dimensional flow, based on 
these assumptions and that 


T, 
H+1= T. (H,+ 1), (9) 


where the suffix “i” refers to a corresponding incompressible case, and 
T,=T.{1+4r(y—1) M’}, 


T, being the recovery temperature and r the recovery factor, which is taken as 0°89, 
and M the external Mach number. 


Alternatively, from a curve given by Young“, we may take 


which fits Young’s curve well for values of M between | and 4. Equation (10) agrees 
with equation (9) if r is taken as 0-8 and H, as 1°5S. 


The skin friction at the surface (7) is given by 


( 


pU? Pe bm (I 


where D and N have their incompressible values, the suffix ““m” refers to a “mean” 
condition relating compressible to incompressible flow, and the suffix “e” refers to 
flow outside the boundary layer. Spence takes 


(T, —-T.)=0:28T. + 0:72T, 
in the case of zero heat transfer, and so 
Tn=T,. (1 +0°128M?*). 


Spence also takes N=5, D=0-0088 and assumes that in the compressible 
momentum equation H, may be taken as constant and equal to 1:5. 


Spence’s’”? equation then becomes 


M®+°2G (M)=0-0106 M*F (M) d(*) +constant, . (12) 


where the suffix 0 refers to an isentropic stagnation point. In this equation it may, 
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TABLE I 


Zero Heat Transfer | Constant Wall Temperature se 


4 1:8 (T,,/T)+22 
7) 
(T,/T,7* (T,/T,)*** 


1+0°128 M? } 


1+0-2M,? 
‘1+0-2M?2 


however, refer to any point in the external flow that can be reached isentropically 
from just outside the boundary layer. a, is the velocity of sound at the reference 
station and B, F and G are found from Table I. 


In this table 7, is the temperature at an isentropic stagnation point. 
If we write M=U./a, and note that 
a./a,=(T./T,)'”, 


Spence’s equation becomes, in the case of zero heat transfer, 


(*) “US? =0-0106 (£ U.S (*) + constant, 


or, expressing U,, 6 and x non-dimensionally in terms of U, and c, 


where R=U,c/%,. 


The axi-symmetric analogy of Cooke asserts that three-dimensional flow along 
a streamline may be calculated as though it were flow about an axially symmetric 
body of radius r, provided that the cross-flow in the boundary layer is small. r is 
connected with p of Section 4, the relation being 


and so, in this case, 
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The momentum equation used by Spence, namely 


d6 du, _ 


becomes, in axi-symmetric flow over a body of radius r, 


dé 
dx 


6 du, 6dr 
° 


The effect of this is to add an extra multiplier r'* in the left-hand side of 
equation (12) and also in the right-hand side inside the integral, while d (x/c) must 
be replaced by d (s/c), since the integral is along a streamline. 


Using equation (4), we thus obtain 


where U., u., 6 and x are expressed non-dimensionally in terms of U, and c. This 
should be compared with equation (5). 


If there is heat transfer and the wall has a constant temperature 7,,, equation 
(13) is altered by the changed values of B, F(M) and G(M). The result is written in 
full in Section 8, namely equations (21), (22) and (23). 


If it may be assumed that there is no leading edge shock, which should be the 
case if the sweep is such that the leading edge is subsonic, then infinity may be taken 
as the reference condition. If, thus, T,, v, are replaced by T., v~. it can be shown 
that, if the pressure coefficient is given by 


Pox 


Cp= 


Te (say), 
_ 

U. _ 
Mx 


M? 5, 


[G+ cos? A 5K)] 


= U,? —U,,’ sin? A= 
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7. Momentum Thickness at Transition 


In equations (5) and (13) the constants are determined by the assumption that 
the streamwise momentum thickness at transition is continuous. Hence it is 
necessary to know this thickness at the end of the laminar part of the flow. As 
already mentioned, the boundary layer over the laminar part may be calculated by 
the independence principle together with some method such as that of Thwaites®”, 
if incompressible, or by Crabtree’s method’ if compressible (see the Appendix). 
In either case a chordwise momentum thickness is calculated fairly simply, but to 
calculate the streamwise momentum thickness is much more difficult. Fortunately 
the former will be a sufficiently good approximation to use in order to start the 
turbulent calculations after transition, if incidence is small, as we shall now show. 


By the use of laminar methods, the value of the “chordwise” momentum 
thickness, namely 


Pelle 


is found, while for the turbulent calculation the “streamwise’”” momentum thickness 


Pelke 
is required. 


If the axes (x, y) (x’, y’) are taken in the same sense (both right-handed or both 
left-handed) and 8 is the angle between the two axes of x (see Fig. 1), we have 


u=ucos sin B . 17) 
v= —usin cos ‘ (18) 
If v’ is assumed to be small in the laminar flow, equation (18) gives 
v utan B, 
and so, substituting in equation (17), 
uo =~ usec B. 


Also, since v,’=0, 
U.=u, sec B, 
an equation which is exactly true. 
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Hence, from equations (15) and (16), 
6 =~ 6, 


and so the chordwise momentum thickness 6, may be used as an approximation to 
the starting value of the streamwise momentum thickness 9, after transition. 


The same arguments show that ~ 8, 


8. Summary of Calculation Methods in Turbulent Flow 
Writing equations (5) and (13) in full we have, in the incompressible case, 


and in the compressible case, for zero heat transfer, 


T. 1:29 1-2 
76 T = Vex 


where the suffix “T”’ refers to transition values, x is measured over the surface of the 
wing normal to the leading edge, and R=U,c/¥,. U., ue, 6 and x are here expressed 
non-dimensionally in terms of U, and c. 


If there is no leading edge shock, 7,, U, and v, may be replaced by T.., U,, and 
Yen and if c, is known, 7./7T,,, M, U./U, and u./U,, may be found from 
equations (14). 


For a wall of constant temperature 7,,, equation (20) is replaced by 


B-1 1:291-2 B-1 1-2 1:2 
U Ue 6 = Uer Oy 


1-244 
Tn 2 7 
R=Uccly. 


344 The Aeronautical Quarterly 


x 


TURBULENT BOUNDARY LAYERS 


The suffix “s” here refers to an isentropic stagnation point in the free stream. 
U., Ue, 9 and x are expressed, as before, non-dimensionally in terms of U, and a 
length c. 


The local streamwise skin friction coefficient c,; may be found from equation (11) 
once @ has been calculated. 


In terms of non-dimensional quantities, as before, 


Tw 
C= 


2-178 


2 


For constant wall temperature 7,,, the expression (1 + 0:128M?)-°** is replaced 
by T,,,/7T. from equation (22). 


For incompressible flow the last two factors of equation (24) are dropped, but a 
more accurate value may be found by means of the Ludwieg and Tillman formula 


—0-678H —0°268 1-732 

0:246 x 10 ‘ 
again expressed non-dimensionally, provided that H has been determined by 
equations (7) and (8). 


The local momentum thickness 5* is found approximately from 
= 


where H is found from equations (7) and (8), or (9) with H,=1-5, r=0-°89, or from 
equation (10). 


9. Discussion 


In the methods given here it is assumed that the independence principle may 
not be used, and so the integration proceeds following a streamline and assuming 
small cross-flow, so that the axisymmetric analogy” may be used. The calculations 
give the streamwise momentum thickness and the “ streamwise displacement thick- 
ness.” If the cross-flow is small, the latter may be taken as a first approximation to 
the actual displacement thickness. It is also to be understood that in deriving the 
equations the skin friction component in the streamwise direction is the same as it 
would be in two-dimensional flow. 


No special equation for H is given in the compressible case, and it is supposed 
that equations (9) or (10) give a sufficiently accurate value. 


In the incompressible case a method for finding H is given, but it seems to be of 
doubtful validity. However, the method gave fair agreement with experiment in one 
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case where the flow was approximately of the sane nature as that considered in this 
paper (see Fig. 3). 


It must again be made clear that this paper applies in the main to flow at small 
incidence and with small cross-flows. *Nowhere must separation be approached. 


Appendix 
LAMINAR COMPRESSIBLE FLOW 
Crabtree’) has pointed out that an approximate solution of the chordwise equation 


may be obtained by assuming the independence principle to hold and using the 
Stewartson-Illingworth transformation as in Ref. 1. 


This leads to the equation 


x 


-3 ( 5 
~6 5 
6,7=0-45y, (=) Uu, T. u,>dx, 


0 


and, if u.=0 at the start (e.g. for a rounded leading edge), 


(=) 


Alternatively we may again use the axisymmetric analogy and integrate along a 
streamline. 


Since in axially symmetric flow, as shown in Ref. 1, 


s 
2—(). -277 -6 27] Sdc 
6?=0-45y, (7) ru, r’USds, 


0 
we have, putting r=Au,/U, and using equation (4), 


-3 
2—().45, -277 -4 4 
6?=0-45y, (7) T, u,U,*dx. 


0 
Expressed non-dimensionally, as before, this equation may be written 
0-45 ( T.)" 
4,, 262 = 4 
U.4u,76 R UAu, dx; 
0 


it has been written in this form so that it may be compared with equation (20). 
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A Small-Bore Pick-Up for the Measurement 
of Fluctuating Pressures 


R. E. FRANKLIN* aad R. B. ARCHBOLDt 


SumMaRY: The paper describes the development of a small-bore probe 
microphone, the frequency response of which is smoothed by arranging that the 
probe tube is terminated in its characteristic impedance. 


1. Introduction 


The work described in this paper originated in a programme of research on jet 
noise at the University of Southampton. Most of the work done at the University of 
Southampton (and indeed elsewhere) on the noise from model jets has been done 
with standard microphones. This is not entirely satisfactory for the following 
reasons. When scaling the results obtained from model jets in an attempt to 
compare them with results obtained from similar experiments on full-size engines, 
the fact that the pressure-sensitive device, namely the microphone, is the same size 
in both cases raises the question of the effect of this on scale factors. Again, the 
results of an experiment in which the diameter of the microphone diaphragm is 
almost half that of the jet are immediately suspect, .or it is not reasonable to assume 
that a microphone of this size does not affect the geometry of the sound field, 
particularly if measurements are carried out quite near to the jet. Ideally the 
measurements taken should be point measurements of the pressures in the sound 
field and there is clearly a need for a very small microphone for these model tests. 
Such a microphone wouid have wide application outside the field of jet noise work; 
for example, it would be useful for the study of pressure fluctuations in boundary 
layers, static pressure fluctuations in turbulent flow, and many problems quite outside 
the field of aeronautics. 


The need for a small diameter microphone was felt most acutely at Southampton 
University when it was decided to initiate experiments on the space-correlation of 
the pressures in the near noise field of a model jet. In these measurements two 
microphones are used, one of which is held stationary at some reference point in the 
pressure field while the other is tracked past it. For certain field separation distances 
the signals from the two microphones are multiplied together and the average value 
of their product is computed. The result of this averaging is called the correlation 
of the two signals and a graph of this correlation against the separation distance is 


*Department of Aeronautics, University of Southampton. Now at the Oxford University 
Engineering Laboratory. 
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called the correlation curve. Examination of similar work” carried out on a full- 
size engine confirmed an estimate made from purely dimensional considerations that 
the correlation was likely to vary rapidly in distances of the order of one jet diameter. 
As a result of this it was clear that, for reasons of physical size and accuracy of 
measurement, the largest microphones which could be used with a 2 in. diameter 
model jet would be of 0-125 in. diameter. An examination of the literature showed 
that the best solution would be obtained by fitting probe attachments to the standard 
microphones and the most promising of the attachments described in the literature 
seemed to be that developed at the Post Office Research Station for use with a 
moving coil microphone®: *?. 


This paper reports a joint investigation carried out at the University of 
Southampton and the Post Office Research Station into the construction of a probe 
attachment for a condenser microphone. 


2. Principles of Probe Attachments 


The problems of design and construction of very small microphones are 
considerable and, even if they are overcome, one is confronted with the difficulties 
associated with very low sensitivity. As an alternative, one thinks of using a small 
bore tube to pick up the pressure fluctuations and then sensing the pressures at the 
far end of the tube with a standard microphone. This is the principle of the probe 
microphone. 


The probe attachment introduces the question of the transmission of sound 
along pipes. In passing down a pipe a sound wave has to overcome the inertia of 
the air in the pipe and its stiffness or acoustical capacitance. Both of these are 
continuous functions of position along the pipe and so the pipe is a system with 
distributed constants. The distributed effective mass, stiffness and resistance to flow 
combine to form a certain impedance, called the characteristic impedance, which is 
the impedance at any point along an infinite pipe of the same dimensions. There 
are no reflections in an infinite pipe; further, if a pipe of finite length is arranged to 
end, or terminate, in another acoustical element, the impedance of which is equal to 
the characteristic impedance of the pipe, then to a sound wave approaching the end 
the pipe is effectively infinite and again no reflections occur. If the pipe is terminated 
in an impedance different from its characteristic impedance, then reflections occur at 
this impedance discontinuity, the reflected wave being positive or negative as the 
terminating impedance is greater or less than the characteristic impedance of the 
pipe. Under these circumstances standing waves can be set up at certain regular 
frequencies; these frequencies are the resonant frequencies of the pipe*. 


It is now clear that a microphone sensing either the pressure in the wall or that 
at the far end of a probe tube will, in general, indicate pressures different from those 
at the open end. This is not in itself a difficulty, since the microphone and probe 
can be calibrated as a unit in terms of the sound pressure at the open end of the tube. 


*A pipe is the acoustical equivalent of an electrical transmission line. 
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The pressure differences can become troublesome in practice, however, if they 
manifest themselves as undulations in an otherwise smooth sensitivity / frequency 
characteristic and this is what usually happens. There are several ways of over- 
coming this. If the tube is made long and its far end is filled with distributed 
absorbing material, the wave is absorbed before any reflections occur and the 
pressure sensed with a microphone in the tube wall will be such that it is simply 
related to the pressure at the open end of the probe, i.e. a smooth characteristic will 
be obtained. Probes built along these lines are described in Refs. 2,4 and 5. They 
are usually of fairly large bore and are not suitable for the applications mentioned 
in Section |. Alternatively, if the microphone is used to sense the pressure at the 
far end of the tube, the resonances which occur may be damped out by placing 
damping material in the tube. This naturally results in a greater loss of sensitivity 
than that due to attenuation in the simple tube. 


If the standard microphone is to be used to sense the pressure at the far end of 
the tube, then one has to provide a cavity coupler; the probe tube and microphone 
are fitted into its opposite faces. The cavity will have a certain impedance 
depending upon its size and, further, the microphone diaphragm will have a certain 
equivalent acoustical impedance. It seems possible that, by altering the shape and 
size of the cavity, the combined impedance of the cavity and microphone might be 
made equal to the characteristic impedance of the probe tube. This method was the 
one followed in Ref. 3 with very satisfactory results. The work reported in Ref. 3 
was mainly concerned with a moving coil microphone; the present work extends this 
to the condenser microphone. 


Ficure 1. Details of 
probe attachment. | 
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3. Construction of the Probe Attachment 


The microphone used at the University of Southampton for standard noise work 
is the Bruel and Kjaer Condenser Microphone Type 4111, and it was proposed to 
build a probe attachment to fit this. The final design is as shown in Fig. 1. 


In this design the brass body of the attachment has been made fairly thick in 
order to cut down the transmission of noise through it. For this reason, too, the 
actual probe tube, which is of stainless steel hypodermic tubing of 0:031 in. inside 
diameter, is enclosed for most of its length by a thick-walled brass tube; the space 
between the two tubes may be packed with insulating material, should tests prove 
this necessary. These precautions are important, in general, because the actual 
pick-up area of the probe tube is so small; they were particularly important with 
reference to the application for which the microphones were intended, namely on 
model jet noise rigs where the ambient noise level is high. 


PROBE TUBE A PRESS FIT 
| 1-416 IN.DIA. 


L 

1-251N. DIA. ——__ 


DRILLED 0-03 
TO MATCH PROBE TUBE 


FIGURE 2. 


A larger scale sketch of the steel coupler is given in Fig. 2. The Bruel and 
Kjaer microphone has a flat diaphragm which stands roughly 0-005 in. proud of the 
capsule. Thus the cavity is cylindrical and of depth d; it is this dimension which is 
altered (by grinding face A) in adjusting the terminating impedance of the probe 
tube. The coupler is pressed on to the face of the microphone by the sponge rubber 
insert shown in Fig. 1, the contact face between coupler and cavity being lightly 
smeareu with Vaseline to prevent air leaks. 


4. Calibration of Probes 


The probes were calibrated on the Rayleigh Disc Resonance Tube at the Post 
Office Research Station. In this apparatus the particle velocity at the middle of a 
resonance tube is measured by the angular deflection of a small, light disc 
suspended there. One end of the tube is closed by a small loudspeaker and the 
other by the microphone whici; is being calibrated. The frequency of the supply to 
the driving loudspeaker is adjusted until simultaneous maxima of the disc deflection 
and microphone signal are obtained, indicating that the tube is in resonance. The 
deflection of the disc then gives the velocity at the velocity anti-node and from this 
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the pressure at the microphone (which is at the pressure anti-node) may be 
calculated. In practice, once an appropriate frequency has been found, the power 
to the loudspeaker is adjusted to give a standard deflection of the disc; this results 
in a standard pressure being applied to the microphone and, from its open circuit 
e.m.f., the sensitivity of the microphone is known at once. 


With this apparatus, the microphone could be calibrated accurately (i.e. to 
within | /10 dB.) over the range 708,000 c./s. 


5. Results of Tests 


The initial cavity depth was fixed arbitrarily at 0-010 in. and the results of the 
calibration with this cavity are shown in Fig. 3, where the resonant characteristics of 
the probe tube may clearly be seen superimposed on a curve which droops because 
of the friction in the tube and cavity. 


The technique was now to reduce the cavity depth by grinding material off 
face A and then to repeat the calibration to determine the effect of the change in 
dimension. From previous work on the moving coil microphone” it was expected 
that, as the cavity depth was reduced, the peak-to-trough ratio of the resonances 
would gradually decrease until a smooth curve was obtained, indicating that the 
probe tube was terminated correctly. In Ref. 3 it is also shown that further decrease 
in the cavity depth causes the peaks to be replaced by troughs and vice versa 
(i.e. the terminating impedance is reduced to less than the characteristic impedance). 


With the condenser microphone it was found that, taking cuts off face A of 
0-001 in. at a time, the changes from one calibration to the next were very small 
until the combined impedance of the microphone diaphragm and cavity closely 
approached the characteristic impedance of the tube. These small changes are 
probably due to the fact that most of the combined impedance is that of the micro- 
phone diaphragm, which is very high, and changes in the cavity impedance have little 
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effect until the critical stage is reached. As this stage was approached, the cuts 
taken from face A were gradually reduced, finally to 0-0001 in. 


The results of the fifth calibration are also plotted in Fig. 3 and illustrate the 
two main effects of the changes in cavity depth, namely that the peak-to-trough ratio 
of the resonances is reduced and that the response curve is raised at the higher 
frequencies. 


The final calibration is shown in Fig. 4. For this calibration the cavity depth is 
estimated at about 0-0015 in., the actual value being difficult to measure because of 
the danger of puncturing the microphone diaphragm. 


With such a drooping characteristic a microphone is unsatisfactory but, if the 
curve of Fig. 4 is examined, it will be seen that the general slope is about —6 dB. 
per octave. This may be overcome by designing an amplifier having a gain 
characteristic rising at +6 dB. per octave to compensate for the fall, and this has 
been done. The response curve of this amplifier is shown in Fig. 4 and the third 
curve is the overall response of the probe microphone and amplifier combination. 
This overall response is almost the same as the frequency response of the original 
capsule and the gain of the compensating amplifier is conveniently such that the 
sensitivity is also unaltered. 


Finally, it should be noted that the smoothed characteristic of the probe 
microphone could be appreciably improved at the higher frequencies if a current 
model of standard microphone were used. It so happened that the model available 
was specifically designed to have a fall in pressure sensitivity above about 2,000 c./s. 
in order that the free field sensitivity could be constant with frequency, thus 
introducing some handicap in the present investigation. With current models 
advantage could be taken of microphones with constant pressure sensitivity 
extending over 10,000 c./s. 
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6. Concluding Remarks 


It has been shown that a probe attachment of the sort described in Refs. 2 and 
3 can be designed for a standard condenser microphone to give a very smooth 
frequency response, which can be readily equalised to give an overall sensitivity 
substantially independent of frequency. 


For the correlation work mentioned earlier it was necessary to produce a 
matched pair of probes. This was successfully carried out, the two microphones 
having responses which differed by less than 0-5 dB. over the useful frequency range. 
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The Aerodynamic Forces Associated with 
Harmonic Motion of Slender Wings on 
Stationary Bodies of Revolution 


R. D. MILNE 


(Queen Mary College, University of London) 


SumMaRY: The slender-body theory is applied to the determination of the 
unsteady aerodynamic forces on slender wing-body combinations when the body 
is stationary and the wing is deforming harmonically. The uniform body is 
circular in cross section, the wings are placed symmetrically and the wing 
deformation mode is necessarily symmetrical about a vertical plane through 
the body axis. The frequency parameter is restricted to that range for which 
the “cross-flow” potential equation is Laplace’s equation. The case of an all- 
moving slender control surface on a body is treated in detail and numerical 
results are given for the forces and moments on body and control surface when 
the control plan form is triangular: the presence of a gap is neglected. 


1. Introduction 


In Ref. | the slender body theory is applied to the calculation of the unsteady 
aerodynamic forces on deforming, slender wings and heaving and pitching, rigid, 
slender wing-body combinations. This paper deals, in an analogous manner, with the 
case of the deforming, slender wing on a fixed body of revolution. The necessary 
conditions on slenderness, frequency parameter and Mach number for which the 
theory is applicable are briefly set down. 


It is an essential restriction of this analysis that the wing motion should be 
symmetrical in the aeroelastic sense : in addition, in the region of the wings the body 
is required to be uniform. 


Thus for a prescribed harmonic wing deformation the generalised forces on the 
wings and the associated lift and moment on the wing-body combination may be 
determined. 


By neglecting the presence of a gap between wing root and body the analysis 
covers the interesting particular case of a combination consisting of the body and a 
slender wing with hinged control surface or all-moving control. This case is 
developed in some detail and, for the case of an all-moving control of triangular 
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plan form, numerical results are given for the lift and moment on both the 
combination and on the control due to a harmonic control rotation. 


The analysis of this paper, with that of Ref. 1, allows a study of the stability of 
the general longitudinal motion of a slender wing-body combination with flexible 
wings and, by an application of the Laplace transform, the calculation of its response 
characteristics. These applications may be made only to motions which are 
sufficiently slow that they do not violate the restrictions placed on the frequency 
parameter. 


NOTATION 
U_ free stream velocity 


p free stream density 
a_ speed of sound in the free stream 
M Mach number of the free stream 
p, free stream pressure 
Oxyz_ right-handed Cartesian co-ordinates 
t time 
@ perturbation velocity potential 
» angular frequency 
v=wc/U, frequency parameter 
c reference length for the combination 
s(x) semi-span of the combination 
p local pressure 


Ap differential pressure across the wing 


X=y +iz =ce”, X -plane 


X,=y,+iz,=§,e", X,-plane 


+ X,-plane 
U,W3;7,,W13V2,W2 perturbation velocity components 
f(x,y) deformation mode of the wing 
R_ body radius 


$, semi-span of the combination at the wing trailing 
edge 


c=R/s, 
r=R/s 
$,=s—R?/s 


b=s,/(2R) 
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Xa, Xi,%_ xX co-ordinates of hinge axis, intersection of wing 
leading edge and body, and wing trailing edge, 
respectively 


Ln4, COefficients in the series expansion, equation (14) 


ze'’' instantaneous wing incidence 
S=s,c, area of the “basic” wing 
A aspect ratio of the wing 


Lw,Mw; lift and moment on the wing, body, and 
combination, respectively 


Q; generalised force associated with the j° wing mode 
f(x,y) j@ wing mode 


L., Ls. Ls 
M.. Ma uc} force and moment derivatives for the combination 


F (c),G(c),H(c),J (+) defining functions for the combination derivatives 


bs li 
Ma, Ma, Ma 


} force and moment derivatives for the wing 


f(r), defining functions for the wing derivatives 


Fn+,(b) see equation (19) 


Note: The use of a bar over a symbol denotes that it is the amplitude of the 
harmonically oscillating quantity represented by the symbol itself. 


Lift is defined as positive upwards, moment as positive nose-up, and incidence 
as positive nose-up. 


2. The Slender-Body Theory 


The configuration under consideration is shown in Fig. | and it is assumed to be 
immersed in a uniform stream of velocity U and density p. Right-handed Cartesian 
axes are used, with the origin O fixed at a point on the body axis and with the x-axis 
parallel to the stream and the z-axis upward. The wing is situated symmetrically on 
the uniform body whose cross section in the region of the wing is circular. The wing 
is represented as a thin plate oscillating about its mean position in the plane z=0 
but always lying in the immediate vicinity of the plane. The trailing edge of the 
wing is straight and conditions downstream of this section are not considered. 
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Ficure 1, Wing-body combination. 


The perturbation velocity potential satisfies the equation 


1 
The conditions holding at the surface of the wing-body combination are that the 
(rigid) body is at zero incidence and the wing is in motion in a specified manner. 


When the wing motion is taken to be simple harmonic, of angular frequency », 
equation (1) becomes 


It is now assumed that the wing is slender, in the sense that 


where s is a characteristic width and c a characteristic length of the configuration. 
It is further assumed that the frequency of the motion is such that the relation 


(Ms/cy’v <1 


is satisfied, where v=wc/U, the frequency parameter, is O (1). 
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Under these restrictions the potential equation (2) reduces to Laplace’s equation 
in two dimensions”, 


The co-ordinate x enters the solution only as a parameter and the boundary condi- 
tions are applied in every “cross-flow plane”, x=constant. The solution in every 
cross-flow plane is complete in itself so long as the total source strength in each 
cross section is zero, such as is the case in antisymmetric flow problems. 


For the configuration considered (having a uniform body) the pressure is given 
by the linear relation 


- =U¢.+% 


or, for harmonic motion, (? = UG,+iwd (5) 


3. Solution for the Velocity Potential 


The cross section in the plane x=constant is that of Fig. 2(a); this physical 


plane will be referred to as the X-plane, where X = y + iz=(e". 


z z 
y¥-PLANE X\-PLANE 
+2R 
-2R 
| 
@) (b) 
FiGurReE 2. 


The boundary conditions on the potential ¢ in the plane x=constant are : — 
(a) ¢, 


(b) On the surface of the body ((=R, O0=4<2r) the radial 
velocity v,=0. 


(c) At any point on the wing (z=0, R<|yv| <s(x)) the prescribed 
normal velocity must be equal to the normal derivative of 4. Let 
the vertical displacement of the point (x,y) on the wing be 
represented by f (x, y) e’*’; then 


w (x, y)=we-“'=¢, (x, y, 0) =iwf (x, y)+ Uf. (xy). (6) 
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(d) Outside the configuration and outside the wake, ¢ (x, y, z) must 
be continuous and, since it is antisymmetric in z, must satisfy the 
condition 


$(x,y,0)=0, |y| > s(x). 
(e) The flow is symmetrical about the plane y=0, that is 
f (x.+y)=f (x, —y), 
$(+y,2)=9(—y, 2) 
and v (0, z)=4¢, (0, z)=0. 


is used to map the cross section of Fig. 2(a) onto the X,-plane of Fig. 2(b). In this 
plane the body is represented by the sides of a vertical slit along the z,-axis and the 
wings by that part of the real axis for which 


0< <s—(R?/s)=s,. 


The relation (7) maps both the inside and the outside of the circle (=R onto the 
X,-plane and in this instance the latter is the required transformation. The 
appropriate inversion is obtained by making the points at infinity correspond; that 
is, the positive sign is selected for the square root, giving 


x=5x, . « 


The points X, = +2iR are branch points. 


The boundary conditions in the X,-plane are, correspondingly : — 
(a) $:,,¢.,—>90 as 0<6,<2r. 


1 
(6) v,=0 for y,=+0, O< |z,| <2R. 
W(x, y,)=w [1+ 


= $2, (x, y,,0), OS Ss, (x). (9) 
{d) $(x,y,.0)=0 for ly,|>s,. 
(e) v,=0 ior y,=0 and for all z, 
2). 
The condition (e) thus contains the condition (b), which may therefore be omitted. 


The restriction that the flow is symmetrical about the y,-plane means that ¢ is 
determined solely by conditions on the wing in the X,-plane. To obtain the solution 
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to this “equivalent wing” problem, it is convenient to introduce the transformation 
X,=s, cosh X,=s, cosh (y,+iz.), . (10) 


in which z, may be taken to vary between 0 and 27 and y, from 0 to 0. The 
X,-plane is thus mapped onto a semi-infinite rectangular region of the X,-plane of 
width 2x. The upper and lower surfaces of the wing are represented by the z,-axis 
between, respectively, 0 and =, * and 27, the correspondence between points on the 
wing being given by 


The remainder of the real axis is given by z,=0, =, 2%, y,=0—> 00. 
The body occupies those portions of the lines z, = 7/2, 37/2 for which 
0< y, < sinh~'(2R/s,). The points at infinity correspond. 


The boundary conditions in the X,-plane are, correspondingly : — 
(a) dv, —>0 as y, > ©. 
(c) 7, (x, 2,)=W, 5, COs Z,) 5, sin Z, = $y, (x, 0, 2,). . il) 


Combining this with the relation (9), 


b cos z,| 
(1 + b* cos? z,)' 


S, Sin 
(12) 


(x, 0, Z,)= Ww (x, 5, COs { 


where b=s,/(2R). 
$(x. y=0, I, 2. 
w,=0 for z,=%/2, 3z/2, 
2y+1)% 2y+1)% 
1 


The potential equation in the X,-plane is 


for which a general solution is 
(K cos Cz, +L sin (M cos Cz, + N sin Cz,), 


where C® is a constant of separation. Boundary condition (a) requires that 
M = N=0, while condition (d) requires that K =0 and that C is any positive integer. 
Condition (e) restricts C to the odd positive integers. Thus the solution of 
equation (13) satisfying the boundary conditions (a), (d) and (e) is 


x 


and the L,,,,, are to be determined from the condition (c). Thus, on the wing, 


~ 
(2n + 1) Sin (2n+ 1) 2, =W,5, Sin Z,. 
0 
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Multiplying both sides by sin (2+ 1) z,, integrating from 0 to = and noting that W is 
symmetrical about =/2, we have 


4s, 
—(2n+ 1) | (a, 5, COS + 


b 
(1+ } sin (2n+1)z, sin Z,dz, 
(15) 


On the wing, (x, 0, Z,)= Lang, Sin (2n+ 1) . (16) 
As required, ¢ is antisymmetric across the wing surface. 


4. Lift and Moment due to an Oscillating Rigid Control 
4.1. GENERAL 

The motion to be considered is that due to the whole or a part of the wing 
pitching harmonically about an axis lying in the plane z=0 and normal to the body. 
Such a motion implies that a variable gap will exist between the wing root and 


body; the presence of this gap is neglected. An assessment of its effect on the lift in 
the steady case is given in Ref. 3. 


Let the line x= x, be the hinge axis; then the function f (x, y) takes the form 
for an all-moving control and 


f (x, y)=0 
=(x,—X) 2, . . . (17d) 


for a flap control, where ze‘«' is the instantaneous angle of incidence. It is convenient 
to continue the analysis for the all-moving control only, since the flap control differs 
from this only in that ¢ = 0 for x < x. 


4.2. THE VELOCITY POTENTIAL 
The coefficients L,,,, are given by 


(2n + 1) Long, = Uz [1 + | x 


/2 

b cos z, } 

2 2, 18 

{1+ sin (2n+1) z, sin z,dz, (18) 


The second term of the integrand is evaluated by replacing sin (2n+1)z, by a sum 
of odd powers of sin z, and using the substitution (1 + 5? cos? z)=u’, which leads to 
integrals of the type 


m 


i (m? du, where m?=b? +1. 


1 
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Writing Lions, = Uz + (x-x)] Fans, 


the first three functions F,, F,, F; are given by 


5[1+ -1)} 


where B=cos~!(1/m). 


1 (ms = 
It may be noted that Lim 5( 5 -1) =2/3 


b+ 
and Lim -1) =%/2. 


The functions F,, F, are zero for R/s=0, 1 and reach their maximum 
numerical values for R/s =~ 0-2. At this value of R/s the relative magnitudes of the 
F.n4,, and thus of the L,,,,, are 


It is seen that convergence is very rapid. 


As a general point, it is of interest to note that, if the L,,,, are regarded as 
Fourier coefficients in an expansion for the function g (z,) sin z, (cf. equation (15)) 
in the range 0 < z, <_=/2, then, since the coefficients in the series for @ are these 
Fourier coefficients divided by 2m+1, the potential is relatively insensitive to the 
harmonics in the Fourier expansion of g (z,) sin z,. 

4.3. Lift AND MOMENT ON THE CONTROL 


The differential pressure across the control is given by 
AB (x, (.+ 5) 
p(x, y)=2pU (6 u? 


and the lift and moment about the origin for both controls by 


Ap dydx 


-M,=2)| xAp dydx,. 
where the integration is over one wing panel. 


November 1960 363 


| | | 


These integrals are most conveniently evaluated in the X,-plane; thus, 


bcos z, 


— sin z, dz,dx, 
(1 + b? cos? z,)'/? } <2 


with a similar expression for M,. 


On substituting for Ap in terms of ¢ in the series form (16), it is seen that the 
lift on a spanwise strip depends on integrals of the same type as those determining 
the L,.,, (equation (18)). As a result, the relative magnitudes of the contributions of 
the successive terms of the series (16) to the lift on a spanwise strip are in the 
proportions 

(0,7 : 


which in numerical terms (for R/s ~ 0-2) are 
100:1/2:1/5. 


It is clear that for practical purposes the lift and moment on the wing may be 
computed using only the first term of the series, when the error is at most 4 per cent. 


4.4. Lift AND MOMENT ON THE WING-BoDy COMBINATION 


The pressure on the body surface is given by 
and the lift and moment about the origin by 


Ly=4| pssin@R dédx . 


B 


—M,=4 | sin 6 R dédx, (24d) 


B 


where the integration is over that part of the body surface for which O0<4< 7/2 
and x, <x < 


Calculation of the lift and moment on the combination is facilitated by the 
following considerations. Let 4’ be the potential (amplitude) due to an isolated 
cylindrical body at incidence 2’e‘~' to the free stream. Then, since ¢ and @’ are both 
solutions of Laplace’s equation (4), the following relation® holds over a closed 
surface S (x) which is taken to be the surface of the body and that part of the plane 
z=0 occupied by the wing in the plane x =constant : — 


= = e 25 
| as dS (x=constant) (25) 


where 7 is the outward normal to S (x). 
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For the potential ¢’, the conditions on S are 
on 


Ons G2 


The downwash w’ ()) is that due to the potential 


=--Uza'sin@, on the body Sz, 


=w’ (y), on the wings Sy. 


, 


= 
evaluated on the plane z=0; thus 
For the potential 4, the conditions on S are those set out in Section 3, namely, 
on 
06 
On ez 


=@ on the body S; 


and =Ww(y) on the wings S,. 


The integral relation (25) taken over body and wings becomes 


2 
| Sin 6 bw dS,. (x=constant) . - 


‘B 


Using this result in the pressure equations (20) and (23), it is clear that the lift 
and moment on the body may be simply related to those on the wing; thus 


Ly=2 | dydx 


. R? 
-M,=2| x Ap-,dydx. . . (285) 


y 


Adding equations (21) and (28), and noting that 


R? 
(1 dyn 


then Lesn=2 | Ap dy,dx=2 | Ap S, sinz,dz,dx . (29a) 


r 
z 


—Mesn=2| x Ap dy,dx=2 | | x aps, sin z, dz,dx. . (29b) 


An immediate consequence of equations (29) is that the lift and moment on the 
wing-body combination depend only on the coefficient L,, so that these results are 
exact within the assumptions of the slender-body theory. The lift and moment on 
the wing or body separately depend on all the coefficients L.,..,. 
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5. The Calculation of Generalised Forces 


The generalised force Q; (for two wing panels) associated with the normalised 
deformation function f; (x, y) is given by 


Apf,dxdy, . . (30) 


where Ap, is the pressure differential across the wing due to the normalised 
deformation function f, (x, y). 


The pressure Ap, is determined from the potential whose coefficients (L.,+,); 
are given by the relations (15) with w (x, s, cos z,)=W, (x, 5, COS Z,). 


The lift and moment on the wing-body combination due to the normalised 
wing deformation f; (x, y) are clearly given by 


| x Ap, dy,dx, . (31d) 


the result (27) again being applied. 


These forces, together with the generalised forces associated with plunging and 
pitching of the rigid wing body combination, allow a general longitudinal stability 
analysis to be undertaken for the combination based on a semi-rigid representation 
of the wing. 


6. Numerical Results for a Triangular Control Surface 
6.1. DEFINITIONS 


In this section numerical results are presented for the lift and moment due to a 
triangular control surface. The origin of co-ordinates is conveniently taken to be 
at the point of intersection of the leading edges when produced to the body axis. 
The triangular wing, so formed by these extensions and by the continuation of the 
trailing edge through the body, will be referred to as the “basic” wing. Let c be the 
root chord and let s, be the semi-span of this “basic” wing: also let e= R/s,. 


Then the root chord of the exposed control surface is 


area of control surface _ 


—— = 2 
and area of “basic” wing il-e7. 


Let x, be the distance of the hinge axis from the origin and x, the distance of the 
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hinge axis from the intersection of the wing leading edge and body; then we define 


* 


so that £,=&+¢ (1-&) 
Let S=s,c be the area of the “basic” wing and let A =(2s,)’/S be its aspect 
ratio; this is also the aspect ratio of the control. 


The notation of flutter derivatives rather than that of stability derivatives 
is used. 


6.2. Lift AND MOMENT ON THE COMBINATION 
Let Lwip=4pU?S {L.+ivLs—v?La} 2 
and Mys3=4pU?Sc {M.+ivM.—v*Ms} 2 


Then F (oe) A 
La=(*/2) F (c)+G A 
Li=(*/2)[H (0) —£,G A 
—M.=(/2) F (¢)—G A 
A 
—Mi=(%/2) (J + (7)] A 


where F(o)=F, Dhan H 


F, (r) dr 


1 
(1 


1 
G =o" | (r) dr, J | 


- F, dr. 


In these expressions r=R/s. In order to carry out the indicated integrations the 
function F, (r) was replaced by the quadratic approximation. 


the maximum error in the integrals is estimated to be 1-5 per cent. The functions 
G, H, J are shown in Fig. 3. 


The movement of the aerodynamic centre of the combination with variation in 
a is shown in Fig. 5, together with the aerodynamic centre position for both wing 
and body at incidence. The function F (c) is the ratio of the lift on the wing-body 
combination when the body is at zero incidence to that when the body is at the 
incidence of the wing: this relation, which is independent of wing plan form, is 
shown in Fig. 5 and agrees with a curve derived by Mirels“ from the complete 
complex potential. 
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Ficure 5. Aerodynamic centre and 
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FiGurRE 6. Aerodynamic centre and lift, control surface. 
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6.3. Lirr AND MOMENT ON THE CONTROL SURFACE 


Let L=4pU*S{l.+ivl,—vls} . (37a) 
and 2. . (37b) 


Then f(c) A 
(c)+ (0) + 


The functions f(c), g(c), h(c), j (7) are given by expressions identical in form to 
equations (35), except that F, (r) is everywhere replaced by (F, (r))’: the functions 
g, h, j are shown in Fig. 4. 


The position of the aerodynamic centre of the control, expressed as a fraction 
of the root chord of the control, is shown in Fig. 6, together with the aerodynamic 
centre position for the control when the body is treated as a reflecting surface. In 
evaluating g (c) in the region 0°8 < z < 1, the approximation (36) is not sufficiently 
accurate for the determination of ¢, and a Taylor expansion of (F, (r))’ about r=1 
was used. 


The control derivatives may be based on the area of the control surface by 
dividing the expressions (38) by (1—<)*. Based on this area, the ratio of the lift on 
the control to that when the body is treated as a reflecting surface is 


f (ce). 
This ratio is plotted in Fig. 6. 


The ratio of the lift carried by the control to that carried by the combination 
is simply F (vc); hence, for small control surfaces, the lift carried by the control 
approaches one half of that carried by the combination. 
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The Reliability of Multiple Rocket Engines* 


GILLIAN RUSTON 


(Rocket Propulsion Establishment, Westcott) 


SuMMaRY: In assessing the merits of multiple versus single engines the possi- 

bility of a loss in reliability due to the use of more working parts must be 

considered. This paper gives a mathematical treatment of this problem and 
numerical examples illustrate the method. 


Introduction 


Suppose that it is required to supply, within a given limited time, a rocket 
engine of given thrust and of acceptable reliability. Will the best solution be a 
single or a multiple engine? Clearly, consideration of weight, efficiency and the 
control system for the missile are among the major issues but this paper is confined 
to a study of the appropriate reliabilities. A simple mathematical treatment of the 
problem is first propounded and then the theory is modified to bring it more nearly 
into line with current knowledge. Each modification is considered separately, and 
the necessary parameters estimated, and then the merits of multiple versus single 
engines are considered on the basis of the suggested modifications taken 
in conjunction. 


NOTATION 

a parameter of thrust-time equation (equation (18)) 

b parameter of development rate equation (equation (19)) 

c parameter of man-power needs equation (equation (21)) 

f (n)=(1 —C’ log. n) 

n number of units in multi-unit engine 

p, probability of failure of single engine giving thrust T 
probability of failure of single engine giving thrust T/n 


probability of failure of single engine giving thrust T/n after development 
time 


probability of failure of single engine giving thrust T after development 


*Correspondence on this paper is especially invited by the Editor. 
Paper received August 1960. 
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q,=1-—p,, reliability of single engine giving thrust T 
q- feliability of combining mechanism 
t development time 
t, actual development time 
t.(n) time needed to combine n sub-engines 
te effective development time 
A parameter of thrust-time equation (equation (18)) 
B_ parameter of thrust-time equation (equation (18)) 
B=T*B 
C parameter of multi-engine combining time (equation (8)) 
C’=C/t 
E number of men who could usefully be employed on a project 
F parameter of man-power needs equation (equation (21)) 


aF 


K 


H_ parameter of development rate equation (equation (19)) 
K number of men available for a project 

L_ parameter of multi-engine combining time (equation (14)) 
P, probability of failure of n-tuple engine giving thrust T 
Q,=1-P,, reliability of n-tuple engine giving thrust T 

T thrust 

a parameter of development rate equation (equation (3)) 

B=a/t 

Nm COefficient measuring effort 

m=K/E, man-power coefficient 


amount of time spent on project _ ; 
= : = priority coefficient 
total working time 


2. Simple Mathematical Model 


Suppose p, is the probability of failure during a run of a rocket engine giving 
thrust T and that p, is the probability of failure during a run of an engine giving 
thrust T7/n. Then the probability of failure of the multiple unit of m of the smaller 
identical engines, total thrust 7, (i.e. the probability of failure of at least one 
engine) is 


assuming that the failures are independent. 
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If np, is small, 


(e.g. if p»=1/1,000, n=4; then np,=4x 10-*, np, —P, <6 10~*), and therefore 
the multiple unit is more or less reliable than the single unit according as 


Suppose that ¢ is the time needed to develop a single rocket engine giving 
thrust T and having probability of failure p (7,1). If T is to be the total thrust of n 
identical engines, the thrust of each is T/m and we may write p(n, t) as the pro- 
bability of failure of one such engine when a time ¢ has been spent on its 
development. A suggested functional form for p(n, t) (to be modified in later 
sections) is 


where 2 and a are constants. We note that 


(i) for fixed m, p(n, t) decreases with increasing t; p(n,0)=1; p(n,t)—>0 
as t—> 0,7 


(ii) for fixed t, p(n, t) decreases with increasing n, i.e. with decreasing engine 
size, 


(iii) for a given probability level ¢ oc (1/n)’, i.e. t OC (thrust)’. 


If the development time f¢ is fixed we have, fom equation (3), 


where 8=2/t is a constant. Using the expression (1), the probability of failure of 
the multiple unit of n engines is 


Taking survival chance as the measure of reliability, 


the reliability of the single engine g,=exp(—§), 


the reliability of the n-multiple unit Q,—exp ( 5 ad 


a 


so that Q,=(q,)"_ 
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Since gq, < 1,Q, 54, according as aS 1; that is, 
if a < 1, the reliability of the multiple unit decreases with increasing n, 


if a> 1, the reliability of the multiple unit increases with increasing n. 


The validity of this conclusion depends, of course, on the validity of equation 
(3) and this may be questioned. Equation (3) may be written more generally in 
terms of thrust as 


where p(T, t) is the probability of failure of an engine giving thrust T after a time f 
has been spent on its development; in later sections modifications to equation (7) 
are considered. 


3. Development of a Composite Engine 


In Section 2 it was assumed that to develop a multiple unit it is necessary simply 
to develop the individual engines making up the composite engine; the relation 


given there implies that the individual engines are completely independent. It has 
been suggested that the simple product rule is unduly pessimistic and that the 
logarithm of reliability should be used in the product rule, instead of the actual 
reliability. This may be accounted for by the human element, perhaps because 
research workers are expected to rise to the challenge of a more complex system, 
but, in the absence of mathematical justification for this supposition, the simple 
product rule is used in this discussion. 


For a multi-engined piloted aircraft it is essential that the engines should be 
completely independent, since failure of one or more engines without loss of the 
aircraft is contemplated and the failure of one engine must not, as far as practicable, 
increase the probability of failure of another. Here this consideration does not hold, 
since failure of one unit means failure of the whole missile to achieve its goal. 
Nevertheless, bearing in mind the need for economy in development, the present 
discussion is limited to a composite engine in which each unit is complete in itself. 


It will still be necessary to construct the n-tuple engine from the single units 
and no allowance was made for this in Section 2, nor for the time needed for 
solving such problems as the synchronising of the firing of the engines. Some of the 
development work for combining units might be carried on simultaneously with the 
development of the individual engines (assuming the necessary man-power to be 
available (Section 6)) but the nature of the work is essentially different and perhaps 
the fairest way of accounting for the development of the combination of units is on 
the time scale. It is suggested that the time needed to develop to an acceptably high 
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level of reliability the composite n-tuple unit from the n smaller units is proportional 
to the logarithm of n, 


t. (n)=C log, n, (8) 


where C is a function of the reliability required. We note that 
() 1.(1)=0, 
(ii) 1, (n) is an increasing function of n, 


(iii) the extra time needed for development in changing from | to 2 engines is 
greater than that needed in changing from 2 to 3 engines and so on. 


The expression (8) assumes that, for instance, the mechanism required for 
8 engines is more than a simple extension of that required for 4 engines, and is 
perhaps pessimistic for the larger multiples. 


If the total time available for development is f, the time left for the development 
of a single sub-engine is 


t—t.(n)=t—C log. n 


and the probability of failure of such an engine is, from equation (3), 


Assuming that the combining mechanism is completely reliable, the probability of 
survival of the composite engines (assuming independence) is 


@ 


where 8=2/t, C’=C/t. We still have the reliability of the single engine giving the 
required thrust as 


q,=¢xp ( —B) 
and so) FE, = (n), 


where f (n)=(1 —C’ log, n) 


Since q, <1, Qn = q, according as f (n)> 1; that is, 
the multiple unit is more reliable than the single unit if f(m) > 1, 


the multiple unit is less reliable than the single unit if f (nm) < 1. 
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Figure 1. Effect on reliability of allowing time for combining sub-engines. 


Values of f(m) are plotted against n for various values of a and C’ in Fig. 1. For 
increased reliability with the multiple engines we must have a > 1. 


It is found, by differentiating (12), that Q, has a maximum value when 


a-—-1-C’ 


(a—1)C’ (13) 


log. n= 
Imposing the condition that a > 1, equation (13) gives a solution for n > 1 only if 


a>1+C’. 


If a<1+C’, then f(n) is a continuously decreasing function for n> 1 and 
therefore in this case a single engine will always give higher reliability than a 
multiple unit. 

To give an idea of the significance of the value of C’, suppose the total time 
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allowed for development is 10 years; then, if C’=0-1, 
t.(2) =0-7 years 
t.(4) years 
t. (32)=3°5 years. 


The main problem of the multi-engine unit will be that of starting the unit and, 
although it may be expected that little in the way of extra equipment liable to 
failure would be needed, it may be thought that the preceding analysis may have 
been unduly optimistic in assuming that the combining mechanism can be made 
absolutely reliable in a limited time. The degree of reliability will depend on the 
time expended on it, and time for development of the combining mechanism may, 
for the purposes of this assessment, be assumed to be obtained at the expense of the 
development of the individual engines. To allow for some degree of unreliability in 
the combining mechanism we may replace (8) by 


L log. n 
log. 


where q. is the reliability of the combining mechanism and L is a constant. The 
behaviour of the expression (14), as m varies while q. is kept constant, is the same as 
that of expression (8) with C replaced by L/log.g., (since 0< q.<1, L will be 
negative). We note that 


(i) t.(0,n)=0, 


(ii) as gq. —> 1, t.(q..n)—> ©, 


t. (q-. nN) = (14) 


(iii) for constant n we have 


t. (Qe, mM) OC log. 


whence qe=exp (= ) 
or p-=1-—exp ( =) 


where A is a positive constant, in agreement with the rate of development given in 
relation (3) for individual engines. 


The probability of failure of a single engine is now, from relation (3). 


log. Ge 


and the survival chance for the composite engine is 


log. 
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It is necessary to find the best division of development time between the individual 
efgines and the combining mechanism by maximising Q, with respect to q.. Purely 
as an illustration, Fig. 2 shows how Q, varies with q., taking as values of the 
relevant parameters 


t=10, n=4, a=1:2, 2=1:0, L= —0:09. 


These values give g, =0-905 and, for g,=0-9, then t,=1-2. 
Differentiating equation (15), 


a 
dq. Qe t 10g, Ge (loge 
and Q, is a maximum when 
(t log. qo—L log. n)? = 


The quantity in the bracket on the left is negative, since t > L log.n/log.q. and logeq. 
is negative. The optimum value of qg. given by equation (16) is plotted against n in 
Fig. 2, taking the same values of the parameters as before in this section. 


OPTIMUM Ge 


0:7 
0-78 0-9 
Figure 2. Variation of the reliability of a 4-unit engine with q,, the reliability of the 


combining mechanism; and of optimum g, with the number of units in a multi- 
unit engine. 
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0:95 
Ficure 3. Varia- 
09 tion of maximum 
reliability of an n- 
fold engine with n. 
0-85 


The maximum value of the reliability is obtained by substituting in 
equation (15) the value of g. obtained from equation (16). Since the reliability of 
the multiple unit is sensitive to the reliability that can be achieved for the combining 
mechanism, two values of the appropriate parameter, L, have been chosen : — 

L=-—0-09 corresponding to g.=0°9 and t,(2)=0°6 

L=-—0009 corresponding to g.=0-99 and (2)=0°6. 
The curves of reliability against number of units are shown in Fig. 3 and, for 
comparison, the broken line is the curve given by equation (10), taking z=1-0, 
C=0°'87, so that ft. (2)=0-6, gq, =0-905 and the combining mechanism is assumed to 


be perfectly reliable. The figure demonstrates the importance of the reliability of 
the combining mechanism if multiple units are to be acceptable. 


4. Me, .udence of Development Time on Size of Projects 


Equation (7) gives the relation 


where ¢ is the development time needed to bring an engine giving thrust T to a given 
level of reliability. The implication of the value of a is suggested by the 
following values 


If it takes 10 years to develop a 200,000 Ib. thrust engine to an acceptable level 
of reliability, the time needed to develop a 50,000 Ib. engine to the same level 
of reliability is, from relation (17), 


2 months if a=3, 24 years if a=1, 
1} years if a=3/2, 5 years if a=4. 
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Ficure 4. Development time against thrust for 1/100 failure rate. 


The value a=1 is the “break even” value of a for which, on the simple theory of 
Section 2, there is neither loss nor gain in reliability for the developt.<nt of a multi- 
unit engine. For a < | the single unit is more reliable than the com”osite engine. 


A comparison of these values with such data on development time as is 
available suggests that the simple power relation 


toc 


is not a fair representation of the facts: there would seem to be a “dead weight” of 
time spent in development irrespective of the size of the project. We therefore 
suggest the relation 


where the values of A and B will depend on the level of reliability required and also, 
to some extent, on when the project is begun, that is, on how much past experience 
is available. Fig. 4 shows curves of time against thrust for three suggested sets of 
values applicable to a given failure rate of 1/100, namely 


A=40 A=30 | A=2°0 
a=0:7 a=0-°55 a=0-43 
B=0-003 B=0-02 B=0-085 


For a project starting at a later date we might expect the values of A and B to be 
less while that of a remains constant. The set 

A=30 

a=0°55 

B=0-02 
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might be modified to 


A=2:0 
a=0°55 
B=0-015 


The resulting curve is shown as a broken line in Fig. 4. 


5. Dependence of Reliability on Development Time 


It seems likely that the rate of development given by equation (7) could be 
improved upon in practice, but that the exponential form is satisfactory. The 
development of any engine will, of course, proceed by steps (if not by leaps and 
bounds) and the functional form can only be expected to represent a mean of 
development of projects of similar size. It is suggested that equation (7) be replaced 
by the expression 


where b is a constant greater than |. The relationship between thrust and develop- 
ment time is then unchanged, but the rate of improvement in reliability is increased. 
Suppose B and a are chosen to relate t and T for a probability of failure of 1/100: 
then H = 0:01 in equation (19). If a certain project has a probability of failure of 
1/100 after 5 years development, Fig. 5 shows the development rate for various 
values of b. 


6. Effect of Limitation of Available Effort 


It would be expected that the time taken to develop a project to a given level of 
reliability would depend on the effort put into the development. This suggests that 
the value of f to be used in equation (7) should be the effective time f. given by 


where #, is the actual time, » is a coefficient measuring the effort put into development 
and0Q<y< 1. 


We further define "=p Nm 


amount of time spent on project 


where »,=priority coefficient = 
y total working time 


number of men working on project 
number of men who could usefully be 
employed on project 


"m= Manpower coefficient = 


and 
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PROJECT HAS FAILURE RATE \/100 
AFTER 5 YEARS —+ 


0-0} 


PROBABILITY OF FAILURE 


‘00! 


20 30 
DEVELOPMENT TIME t (YEARS) 
Ficure 5. Increase in reliability with development time. 


Let E be the number of men who could usefully be employed on the project; 
E may be expected to be a slowly increasing function of the thrust that the engine is 
designed to produce, say 


E(=F(-) , 


where F and c are constants. For example, if the total thrust to be provided by 
n engines is 200,000 Ib. we may have F=100, c=0-3; the resulting curve for E is 
shown in Fig. 6. In practice ¢,, », and 7m E, number of men available (= K), will be 
fixed, so that 


te=Np"mla = Np t, where E(n)>K 


=f, where E(n)=< K. (22) 
In Fig. 6, ¢, is plotted against n for parameter values K = 56, »,=0°9, t, = 20. 


In the region for which E (n) > K, we have, from equation (3), 
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oF 
tanyK 
and when E(n)<K 


where G 


p(n, t)=1—exp 


Comparison of equations (23) and (24) with equation (3) shows that the placing 
of some limit on available effort has, in effect, increased the value of z. This 
indicates a decrease in reliability of an engine of any size at any stage of its develop- 
ment. However, the decision as to whether the single or the multiple engine is 
better depends on the value of the parameter a (Section 2). 


Equation (23) shows that, when less men are available than could usefully be 
employed, the power to which n is raised is increased (in effect, a is increased). We 
now have the criterion: —the multiple unit is more or less reliable than the single 


engine according as a+c=1, which should be compared with that given just after 
equation (6). 


7. Modified Mathematical Theory 


Each of the modifications introduced in the preceding sections requires a change 
in the form of the relation (7) and the estimation of the values of one or more 
parameters—it seems that the simple theory of Section 2 with its two parameters has 
little but its elegance to commend it. ; 


Omitting the question of development effort (Section 6) which introduces a 
discontinuity very special to the conditions of the particular project, we may combine 


100 20 


10 


Ficure 6. Illustration of effect of limitation of effort. 
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Ficure 7. Contours of reliability for a single engine of thrust T after development time ¢. 


the modifications of Sections 4 and 5 and replace equation (7) by 
e. bd 
p(T.)=1-exp { -H(4 


Taking a representative set of parameters 


H=001, A=2:0, a=055, B=0-015, b=3-0, 


contours of p(T, 1) against T and ¢ have been calculated and are shown in Fig. 7. 


The final consideration is the modification of Section 3: if it is assumed that 
the combining mechanism can be made absolutely reliable, equation (25) can be 
replaced by 


A+ BT° (26) 


p(T, t) =1-exp{ F log. n 
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or, setting thrust in terms of n, 


+ An‘ +B’ 
p(n,t)=1 exp { 


where B’=T°B and T is the total thrust of n engines. Then 


An’ +B’ ] 
(t—C log. n) 


Q,=exp { 


and q =exp { } 
Clearly the behaviour of Q, as a function of m depends on the values assumed for 
the various parameters. Fig. 6 shows P,=1-—Q, plotted against n for the sets of 
values given in Table I. 


TABLE I 


Total t 


200,000 OO! 20 0-55 0-015 9-348 10 0-13 “3 30 
200,000 0-01 2:0 0-55 0-015 9-348 20 0-13 3-0 
200,000 001 20 0-55 0-015 9-348 10 0-087 “87 3-0 
200,000 0-01 2-0 0-55 0-015 9-348 10 0-13 2-0 
200,000 0-01 30 0-7 0-002 10°:274 10 0-13 4 30 


If allowance is made for some degree of unreliability in the combining 
mechanism, equation (28) must be modified to give 


An‘ + B’ 
(1. 
log. 


Q.=q.exp J) -H 


Differentiation of this equation shows that Q, will take its maximum value with 
respect to variation in gq. when 


(t log. g.—L log. n)’**  —bHL(An* + B’)’ log. n 
(log. = . (30) 


(which may be compared with equation (16), where 6 was unity). The optimum 
value of q., obtained from equation (30), is substituted in equation (29) to give 
Curves 6 and 7 in Fig. 8 with the sets of parameters given in Table II. 


TABLE Il 


Total t 
H 
6 200,000 0-01 0-55 0-015 10 0-09 30 


7 200,000 0-0! 2: 2-0 0-015 10 0-009 3-0 
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Ficure 8. Probability of failure of composite engine of n units. 


These two curves should be compared with Curve 3 in Fig. 8, for which the 
combining mechanism was assumed to be perfectly reliable. 


8. Conclusions 

It is clear from equation (2) that for a composite engine to reach a given level of 
reliability the sub-engines will each individually have to be more reliable than 
would a single engine designed to give the required thrust. When only a limited 
time is available for development, however, the composite engine may have some 
advantage. 

While these conclusions may seem too obvious to justify the mathematical 
structure presented, it is considered that such a method should enable useful 
predictions on the reliability of future rocket-propelled missiles to be made from an 
analysis of the growth of reliability of past projects. 
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The Quasi-Cylinder of Specified Thickness 
and Shell Loading in Supersonic Flow 


H. PORTNOY 


(Department of Mechanical Engineering, Royal Technical College, Salford) 


Summary: The methods of the operational calculus are used to obtain a linear 
approximation to the shape of the mean camber surface of a quasi-cylinder in a 
supersonic flow in terms of its shell thickness and loading distributions. The 
analysis deals with a generalised quasi-cylinder; that is one which, although lying 
close to a mean cylinder, need not possess axial symmetry. The quasi-cylinder 
is also permitted to be within the small disturbance field of other separate 
components, e.g. a centre-body. Because the linearised theory is inadmissable 
for internal duct flows close to and beyond the first reflected characteristic cone, 
the present solution is likewise invalid close to and beyond the position where 
this characteristic meets the mean cylinder. The work given here enables the 
camber shapes of “ring-wings”, which have been used theoretically to reduce or 
even nullify the wave-drag of a central slender-body, to be found. An example 
illustrates the general method. 


Introduction 


The theory of quasi-cylinders is due to Lighthill®’ and Ward’. Both dealt with 
the external and internal flow problems for axisymmetric quasi-cylinders. It should 
be noted that Ward’s” internal flow analysis forms a correction to that in Ref. 1. 
Ward also considered small-incidence effects. 


The concept of a quasi-cylinder has been extended by Nielsen and Pitts to 
include non-axisymmetric shapes lying close to a mean cylinder and they have dealt 
with the external flow past such bodies. 


In this paper a generalised quasi-cylinder of this kind is considered. In general 
it can be subject to the perturbation field of other separate components together 
with the main stream, which is parallel to its mean cylinder axis in the undisturbed 
state. Its own perturbation field must be such as to satisfy the linearised equation of 
supersonic flow, the conditions at infinity and on the Mach surface from the intake 
edge and, when combined with the other two component flows mentioned, it must 
produce a resultant flow tangential to the outer and inner surfaces of the 
quasi-cylinder. 
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This paper deals with the problem of finding the shape of the mean camber 
surface of a quasi-cylinder of this kind when the thickness and loading distributions 
of the shell are specified, together with the variation of the intake opening radius with 
angular position. The mean camber surface lies midway between the outer and 
inner surfaces. Apart from the usual conditions to be satisfied in quasi-cylinder 
theory®: *, the flow must remain supersonic inside the duct for as far downstream 
of the inlet as the solution is used. Ward has observed thai ile itterna! flow results 
of linearised theory are of doubtful validity as the first reflected characteristic cone 
(which meets the mean cylinder at a distance 2Ba from the inlet) is approached, even 
if the flow is still supersonic. Thus the present work, which assumes the linearised 
supersonic solution to apply within the duct, is valid, at most, for a distance some- 
what less than 2Ba from the inlet. 


The relationships required are obtained, in operational form, from the standard 
solutions of the transformed linearised equation”. On inversion, the thickness and 
loading distributions are seen to enter into the solution through separate superposi- 
tion integrals and with different influence functions, each of which can be expressed 
in terms of complete elliptic integrals. Their effects are therefore simply additive. 


An example of the use of the method is given in Section 4. The mean camber 
surface shape is found for an axisymmetric duct with a conical centre-body, the duct 
having constant shell loading per unit area and linearly increasing shell thickness 
from the inlet onwards. Loading and thickness are interrelated by making the inner 
surface at the inlet lip have zero streamwise slope. Thus the loading at the forward 
edge (and hence everywhere) can be readily found, because two-dimensional theory 
is applicable at the lip. 


The work given here finds, perhaps, its most interesting use in the design of 
“ring-wing”-body combinations of small or even zero theoretical wave-drag. In 
Ref. 5, the method of designing the correct combination of central body of revolution 
with shell thickness and loading distributions of the ring-wing surrounding it is 
considered. The present method enables the mean camber surface shape of the 
ring-wing corresponding to the required distributions of thickness and loading to 
be found. 


NOTATION 
a_ radius of mean cylinder 


a,(z) Fourier coefficient of R (6, z)— R (6, 0) 
A,(p) arbitrary function of p in operational solution for ¢ outside the cylinder 
b, (z) Fourier coefficient of R (@, z)— R (6, 0) 
B= /(M?-1) 
B,(p) arbitrary function of p in operational solution for 4 outside the cylinder 
c.(z) Fourier coefficient of v (6, z)/U 
C,(p) arbitrary function of p in operational solution for 4 inside the cylinder 
d,(z) Fourier coefficient of » (@, z)/U 
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D,(p) arbitrary function of p in operational solution for 4 inside the cylinder 
e, (z) Fourier coefficient of t (4, z) 
E complete elliptic integral of the second kind 
f,(z) Fourier coefficient of z) 
F(z) any function of z 
F(p) the Laplace transform of F (z) 
2,(2),4,(z) Fourier coefficients of / (4, z) 
H(z) the Heaviside unit function. H (z)=1 for z > 0 and zero for z <0 
I, the modified Bessel function of the first kind of order n 
K complete elliptic integral of the first kind 
K,_ the modified Bessel function of the second kind of order n 
1(4,z) radially outward load per unit area of the shell divided by 4U? 
M Mach number 
p the Laplace transform variable 
q.(z) function whose Laplace transform is p /, (p) K, (p) 
Q.(z) function whose Laplace transform is 4p /,’ (p) K,’ (p) 


R(6,z) radius of mean camber surface 


R,.(z) function whose Laplace transform is —4p ipl (p) K, (p)] 


u perturbation velocity vector 
U free stream velocity magnitude 
7 (@,z) radial velocity produced by other components at mean cylinder 
y=2/(Ba) 
a2 semi-cone angle of centre-body in example 
€ constant value of dt/dz in example 
velocity potential such that u=UVo 
p free stream density 
V gradient operator 


2. Operational Solution for the Mean Camber Surface Shape 


The quasi-cylinder is immersed in a stream of initial velocity U, Mach number 
M and density p. Cylindrical polar co-ordinates (r, 6, z) are taken, with the axis Oz 
lying along the mean cylinder axis parallel to the undisturbed stream and the origin 
at the centre of its nose section. The mean cylinder radius is a and the symbol B 
denotes /(M?—1). The sharp edged intake opening has a radius R (4, 0), the mean 
camber surface radius being R (4, z), where 


R (6, z)= [a, (z) cos n6 + b, (z) sin n9]+R(6,0). 
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Other components besides the quasi-cylinder may produce a radial velocity 
distribution at the mean cylinder given by 


v(6,z)=U (z) cos + d, (z) sin n6]. . ; (2) 


The specified thickness of the shell is 


(6, z)= [ea (z) cos +f. (2) sin n6], | 


while the radially outward load per unit area of the shell, divided by 4pU’, is 


1(6, z)= [g. (z) cos n6+h,(z)sinné]. . (4) 


= 


If ¢ is the perturbation velocity potential of the quasi-cylinder, so that its 
perturbation velocity vector u is given by 


then the linear approximation to the pressure coefficient produced by the cylinder’s 
field is —209/0z. Since the pressure due to the other components is continuous 
across the quasi-cylinder, it follows that 


1(6, z)=2 - | 


to a first approximation. 


¢ satisfies the linearised equation’, vanishes far away from the cylinder and 
also on the wave fronts emanating from the inlet edge, and satisfies the tangency 
conditions at the duct surfaces. 


The one-sided Laplace transform is now introduced ©. 


The relationship between a function F (z) and its transform F (p) is indicated 
by the symbolism of van der Pol and Bremmer“? thus : — 


In Chapter 8 of Ref. 4, it is shown that the operational solution of the linearised 
equation suitable for the region outside a quasi-cylinder is 


[A, (p) cos n6 + B,(p)sinn6]K.(Bpr).. . . (8) 


where A,, B, are arbitrary functions of p, and K, is the modified Bessel function 
of the second kind of order n. 


The type of operational solution suitable for the internal region is shown to be 


= = [C.(p) cos . . 9) 
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C,, and D, being again arbitrary functions, while /, is the nodified Bessel'”’ function 
of the first kind of order n. 


The boundary conditions of tangency at the outer and inner surfaces are, 
respectively, to first order, 


Co v(6,z) OR et 
v(6,z) OR Cr 


(10) 
and 


If the operational forms of these are taken and the appropriate series (1), (2), (3), 
(8) and (9) are used, equation of corresponding Fourier coefficents yields 


Bp A, (p) K,’ (Bpa)+C, (p)=p (p)+ 4p (p) , (12) 
Bp B,(p) K,’ (Bpa)+d,(p)=p (13) 
Bp C, (p) 1,’ (Bpa)+C, (p)= p Ga (p) — 4p (p) (14) 
and Bp D, (p) 1,’ (Bpa) +d, (p)=pb, (p)—4p f. (p) 


Also (4) and (6), with (8) and (9), yield similarly 
(p)=2p [A, (p) K, (Bpa) (p) 1, (Bpa)] (16) 
and =A, (p) =2p [B,(p) . (17) 


If A,, B,, C, and D, are now eliminated from (16) and (17) by substitution from 
equations (12) to (15), the following relationships result : — 


4Bg, (p\l,’ (Bpa) (Bpa)=[pa, (p) —C, (p)) (Bpa) (Bpa) — 1, (Bpa) K,’ (Bpa)) 
+4 pe, (p) (Bpa) K, (Bpa)+1,(Bpa) K,’(Bpa)}] (18) 

and 

4Bh, (p) 1,’ (Bpa) K,’ (Bpa)=[pb, (p) —d, (Bpa) K, (Bpa) — 1, (Bpa) K,’ (Bpa)] 
+ 4pf. (p) (Bpa) K, (Bpa)+1, (Bpa) K,’ (Bpa)]. (19) 


If use is made of the Wronskian determinant'’’’ for /, and K, and the formula 
for derivative of a product, these simplify to 


a, (p)= C, (p) + 4B?ag, (p) 1, (Bpa) K.’ (Bpa) — 4pe, [/. (Bpa) K, (Bpa)} 
(20) 


and 
= d. (p) + 4B?ah, (p) I,’ (Bpa) Ky’ (Bpa) 4p}. (p) (1. (Bpa) K,, (Bpa)}. 


(21) 
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equations (20) and (21) may be formally inverted, using certain standard operational 
rules’, to yield 


a, (2) = (u) du+ B (w) du+ 24) 


and b, (z)= 0. (u) du + wo du (25) 


0 0 


These are the required relationships determining the slope of the mean camber 
surface from equation (1). It remains to perform the actual inversions of Q, (p) and 
R, (p) in (22) and (23). 


3. The Functions Q, (z) and R, (z) 
If a subsidiary function q, (z) is defined by 


Gn (Z)= pl. (p)K,(p), . (26) 
then it follows from (23) directly, using a standard operational rule’, that 
To invert Q, (p) we observe that 
tp (p) (p)]=4 (p) Kx” (p) + 21,’ (p) Ka’ (p) + (p) Kx (p)] 
(p+ (1+ %) Kato) + (P) Ke 
+ 4pK, (p) ; I,’ (p)+ (1+ 1.(p)) | 
=-} (p) K, | +4 (1 + Pla (p) K, (p)+ 
where the differential equation for the modified Bessel functions has been used. 
Rearranging and inverting the terms :— 


Q, (z)=($2? — 4) Jove (29) 
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Now, by a formula given in Watson (Ref. 7, p. 441), 
9 
pl,(p) K.(p)= | p K, (2p cos 9) cos 2n6 dé 


and so, by inversion under the integral sign'*’, 


te 
2(-y ~1)cos 


(2)= — j 
z 
20086 -1) 


where H is the Heaviside unit function. Thus 


d6, . 


(-Y cos 2né 


Gn (2) = (32) 


H (1 —j}z)cos— 


It is readily shown from (32) that gq, (z) is expressible in finite series of elliptic 
integrals with modulus 4z for z < 2 and modulus 2/z for z> 2. In particular 


qo (2)= K (42) H(2—2)+ 2 (33) 


K being the complete elliptic integral of the first kind with modulus as indicated. 
Hence, from (29) and (27), 


9 
(E being the complete integral of the second kind) and 


R, (z)= 2K Q-2)+ 1K (35) 

In practice, owing to the limitations on the theory mentioned in Section 1, the 
functions are only needed for z < 2. 

In general it can be shown readily, from (32), that q, (z) has one singularity, at 
z=2, of the form {(-—)**'/(2=)} log |2—z! and so, from (27), R,(z) has precisely 
the same form of singularity at z=2 only and, from (29), Q,(z) has a single 
singularity of the form {(—)"*'/(4=)} log |2—z!. 

Series for Q, and R, for small z and for very large z may be established easily 
by standard operational methods“: *. 


4. Example 


The method is now used to find the camber surface shape of an axisymmetric 
quasi-cylinder with a conical centre-body whose semi-cone angle is 2 and which 
starts at z= — Ba, so that the Mach wave from its nose just meets the intake lip 


(see Fig. 1). 
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Ficure 1. 


The inlet radius R (0) is taken as the radius of the mean cylinder and is denoted 
as usual by a; the thickness of the shell ¢ (z) is given by 


where ¢ is a small positive constant and the shell loading is constant, being fixed by 
making the internal duct surface have zero streamwise slope at the inlet. Since, close 
to the inlet, the flow is two-dimensional it follows that 


2¢ 
B° 


It is readily shown™ that at the mean cylinder the radial velocity induced by 
the central cone for z > 0 is 


As the system is axisymmetric, equation (1) gives R (z)— R (0)=a, (z) and so, 
from (24), (36), (37) and (38), 


0 


(39) 


I(z)=- (37) 


For the present simple case the superposition integrals can be evaluated 
analytically, as can the first integral too. 
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QUASI-CYLINDER IN SUPERSONIC FLOW 


The result of performing these integrations when z < 2Ba is 


{ 1+ [w+ DV 1)? (y+ 1) | + E Gy) - 


40) 


where y has been written for z/(Ba). 


This has been shown in Fig. 1 for «B=0-12B=0-01, as a graph of {R(z)—a}/a 


against z/(Ba). The slope of the line has a logarithmic singularity when z/(Ba)=2 
(which is barely discernible with the scale used). This corresponds, of course, to 
the breakdown of the theory mentioned in Section 1. 
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REVIEW 


Review 


Progress in Solid Mechanics, Volume I. 1. N. Sneddon and R. Hill (Editors). North 
Holland Publishing Co., Amsterdem, 1960. 448 pp. Diagrams. £5. 


The time has long passed since the classical theory of elasticity formed an adequate 
basis for the solution of all of the practical problems of solid mechanics. However, 
knowledge of the progress that has taken place in recent years in the fields of plasticity, 
large deflection elasticity, viscoelasticity and thermoelasticity is not so widespread as it 
ought to be. The series of books of which this is the first volume is therefore very 
welcome indeed. 


This book contains eight short monographs of the highest standard and of great 
scientific interest. The titles are as follows:— I—Viscoelastic Waves by S. C. Hunter. 
II—Matrices of Transmission in Beam Problems by K. Marguerre. II1I]—Dynamic 
Expansion of Spherical Cavities in Metals by H. G. Hopkins. IV—General Theorems 
for Elastic-Plastic Solids by W. T. Koiter. V—Dispersion Relations for Elastic Waves in 
Bars by W. A. Green. VI—Thermoelasticity: The Dynamical Theory by P. Chadwick. 
Vil—Continuous Distributions of Dislocations by B. A. Bilby. VIII—Asymmetric 
Problems of the Theory of Elasticity for a Semi-Infinite Solid and a Thick Plate 
by R. Muki. 


It is not, of course, possible to do justice to this wealth of material in a short review, 
but a few brief indications of content may stimulate interest in particular topics. The 
material basis of Chapter I is expressed by a stress-strain equation in which the modulus 
is complex and frequency-dependent. This is applied to uni-axial stress pulses, compared 
with experiznent and generalised to three dimensions. Chapter II presents a systematic 
development of the Holzer-Myklestad method for the analysis of the vibration of (say) a 
continuous bean) composed of segments of differing elasticity, with concentrated masses 
and supports of various kinds. Chapter III surveys existing knowledge concerning the 
explosive expansion of cavities and the resulting elastic-plastic waves. This difficult 
subject is treated with great mathematical elegance. Chapter IV gives a very thorough 
account of the whole field of small deflection plasticity theory and of limit design. It will 
be invaluable to anyone wishing to do work in this field. Chapter V deals with the 
important and difficult problem of wave propagation in bars and throws much light on 
the adequacy of approximate caiculation and on the way in which the basic elastic waves 
combine to give signal velocities of different magnitude. Chapter VI discusses the basic 
principles of thermoelasticity, following Biot, and applies the results to wave propaga- 
tion in elastic solids, determining quasi-elastic and quasi-thermal waves. Readers should 
notice the misprint in equation (2.1). Chapter VII contains a profound study of 
continuously dislocated crystals and related non-Riemannian spaces, which the author 
believes will lead to a theory of macroscopic plasticity. Finally Chapter VIII deals with 
the solution of elastic problems for plates and half-spaces.—w. s. HEMP. 
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